Anomalous scaling of a passive scalar in the presence of strong anisotropy 
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Field theoretic renormalization group and the operator product expansion are applied to a model 
of a passive scalar quantity 9(t,-x.), advected by the Gaussian strongly anisotropic velocity field 
with the covariance oc 5{t — t')|x — x'|^. Inertial-range anomalous scaling behavior is established, 
and explicit asymptotic expressions for the structure functions S'„(r) = ([S(i,x + r) — ^(t,x)]") 
are obtained; they are represented by superpositions of power laws with nonuniversal (dependent 
(*~^ I on the anisotropy parameters) anomalous exponents, calculated to the first order in e in any space 

• dimension d. In the limit of vanishing anisotropy, the exponents are associated with tensor composite 

' operators built of the scalar gradients, and exhibit a kind of hierarchy related to the degree of 

, anisotropy: the less is the rank, the less is the dimension and, consequently, the more important 

is the contribution to the inertial-range behavior. The leading terms of the even (odd) structure 
' functions are given by the scalar (vector) operators. For the finite anisotropy, the exponents cannot 

be associated with individual operators (which are essentially "mixed" in renormalization), but the 
aforementioned hierarchy survives for all the cases studied. The second-order structure function 5*2 is 
studied in more detail using the renormalization group and zero-mode techniques; the corresponding 
fT^ . exponents and amplitudes are calculated within the perturbation theories in e, 1/d, and in the 

anisotropy parameters. If the anisotropy of the velocity is strong enough, the skewness factor 
83/82^'^ increases going down towards to the depth of the inertial range; the higher-order odd ratios 
increase even if the anisotropy is weak. 



PACS number(s): 47.27.-i, 47.10. -fg, 05.10.Cc 



I. INTRODUCTION 



Recently, much effort has been invested to understand inertial-range anomalous scaling of a passive scalar. Both 
the natural and numerical experiments suggest that the violation of the classical Kolmogorov-Obukhov theory 
i is even more strongly pronounced for a passively advected scalar field than for the velocity field itself; see, e.g., 
[^,|| and literature cited therein. At the same time, the problem of passive advection appears to be easier tractable 
theoretically. The most remarkable progress has been achieved for the so-called Kraichnan's rapid-change model 
' for the first time, the anomalous exponents have been calculated on the basis of a microscopic model and within 
i regular expansions in formal small parameters; see, e^, J^p^ and references therein. 
^ ■ Within the "zero- mode approach," developed in |q-|lO|| , nontrivial anomalous exponents are related to the zero modes 
• ^ ' (unforced solutions) of the closed exact equations satisfied by the equal-time correlations. Within the approach based 
513 ■ on the field theoretic renormalization group (RG) and operator product expansion (OPE), the anomalous scaling 
^ ' emerges as a consequence of the existence in the model of composite operators with negative critical dimensions, 
which determine the anomalous exponents pd)-|ir 



Another important question recently addressed is the effects of large-scale anisotropy on inertial-range statistics 
^ , of passively advected fields [p| JlO|Jl^ -|l9t and the velocity field itself ||2^,|l| . According to the classical Kolmogorov- 
- - ' Obukhov theory, the anisotropy introduced at large scales by the forcing (boundary conditions, geometry of an 
obstacle etc.) dies out when the energy is transferred down to the smaller scales owing to the cascade mechanism 
10,^. A number of recent works confirms this picture for the even correlation functions, thus giving some quantitative 
support to the aforementioned hypothesis on the restored local isotropy of the inertial-range turbulence for the velocity 
and passive fields ||r^ , p^ -|2l[| . More precisely, the exponents describing the inertial-range scaling exhibit universality 
and hierarchy related to the degree of anisotropy, and the leading contribution to an even function is given by the 



exponent from the isotropic shell ]14| , |l5| , |l7| - |l9|21 |. Nevertheless, the anisotropy survives in the inertial range and 
reveals itself in odd correlation functions, in disagreement with what was expected on the basis of the cascade ideas. 
The so-called skewness factor decreases down the scales much slower than expected [pPj9|,|l0t , while the higher-order 
odd dimensionless ratios (hyperskewness ets) increase, thus signalling of persistent small-scale anisotropy p^ , p^ , p^ . 
The effect seems rather universal, being observed for the scalar and vector fields, advected by the Gaussian 
rapid-change velocity, and for the scalar advected by the two-dimensional Navier-Stokes velocity field |p^ . 

In the present paper, we study the anomalous scaling behavior of a passive scalar advected by the time-decorrclated 



strongly anisotropic Gaussian velocity field. In contradistinction with the studies of [l4|,pyiC|Jl4[|, where the velocity was 
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isotropic and the large-scale anisotropy was introduced by the imposed linear mean gradient, the uniaxial anisotropy 
in our model persists for all scales, leading to nonuniversality of the anomalous exponents through their dependence 
on the anisotropy parameters. 
The aim of our paper is twofold. 

Firstly, we obtain explicit inertial-range expressions for the structure functions and correlation functions of the scalar 
gradients and calculate the corresponding anomalous exponents to the first order of the e expansion. We show that 
the exponents in our model become nonuniversal through the dependence on the parameters describing the anisotropy 
of the velocity field. Owing to the anisotropy of the velocity statistics, the composite operators of different ranks mix 
strongly in renormalization, and the corresponding anomalous exponents are given by the eigenvalues of the matrices 
which are neither diagonal nor triangular (in contrast with the case of large-scale anisotropy). In the language of the 
zero-mode technique this means that the SO{d) decompositions of the correlation functions (employed, e.g., in Refs. 
p9|) do not lead to the diagonalization of the differential operators in the corresponding exact equations. 

Nevertheless, the hierarchy obeyed by the exponents for the vanishing anisotropy survives when the anisotropy 
becomes strong. This fact (along with the nonuniversality) is the main qualitative result of the paper. 

Another scope of the paper is to give the detailed account of the RG approach and to compare the results obtained 
within both the RG and zero-mode techniques on the example of the second-order structure function. 

The paper is organized as follows. 

In Sec. OjWe give the precise formulation of the model and outline briefly the general strategy of the RG approach. 
In Sec. pi| , we give the field theoretic formulation of the model and derive the exact equations for the response 
function and pair correlator. 

In Sec. IV, we perform the ultraviolet (UV) renormalization of the model and derive the corresponding RG equations 
with exactly known RG functions (/3 functions and anomalous dimensions). These equations possess an infrared (IR) 
stable fixed point, which establishes the existence of IR scaling with exactly known critical dimensions of the basic 
fields and parameters of the model. 

In Sec. 0, we discuss the solution of the RG equations for various correlation functions and determine the depen- 
dence of the latter on the UV scale. 

In Sec. 0, we discuss the renormalization of various composite operators. In particular, we derive the one-loop 
result for the critical dimensions of the tensor operators, constructed of the scalar gradients. 

In Sec. |VI1| , we show that these dimensions play the part of the anomalous exponents, and present explicit inertial- 
range expressions for the structure functions and correlations of the scalar gradients. 

In Sec. Vin, we give the detailed comparison of the RG and zero-mode techniques on the example of the second- 
order structure function, compare the representation obtained using RG and OPE with the ordinary perturbation 
theory, and calculate the amplitude factors in the corresponding power laws. 

The results obtained are reviewed in Sec. IX. 



II. DEFINITION OF THE MODEL. ANOMALOUS SCALING AND "DANGEROUS" COMPOSITE 

OPERATORS. 

The advection of a passive scalar field 9{x) = 9{t,x.) in the rapid-change model is described by the stochastic 
equation 

VtO^iyoM + f, Vt = dt+v,d,, (1) 

where dt = d/dt, di = d/dxi, vq is the molecular diffusivity coefficient, A is the Laplace operator, v(a;) is the 
transverse (owing to the incompressibility) velocity field, and f = f{x) is a Gaussian scalar noise with zero mean and 
correlator 

{f{x)f{x'))=5{t-t')C{r/£), r^|x-x'|. (2) 

Here I is an integral scale related to the scalar noise and C{r/t} is a function finite as £ ^ oo. With no loss of 
generality, we take C(0) = 1 in what follows. The velocity v(x) obeys a Gaussian distribution with zero mean and 
correlator 

{v,{x)v,{x')) = Do ^i^-^ J dkT,,{k) (fc2 + m2)-'^/2--/2 exp[ik(x - x')]. (3) 

In the isotropic case, the tensor quantity Tij(k) in (||) is taken to be the ordinary transverse projector, Ty(k) — 
Pij(k) = Sij—kikj/k'^, k = |k|, Dq > is an amplitude factor, 1/m is another integral scale, and d is the dimensionality 
of the x space; < e < 2 is a parameter with the real (Kolmogorov) value e — 4/3. The relations 
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Dq/vq EE go EE 



(4) 



define the coupling constant (i.e., the formal expansion parameter in the ordinary perturbation theory) and the 
characteristic UV momentum scale A. In what follows, we shall not distinguish the two IR scales, setting to ~ l/l. 
The issue of interest is, in particular, the behavior of the equal-time structure functions 

SN{r)^{[e{t,^)-e{t,^)f) (5) 

in the inertial range, specified by the inequalities l/A<<r<<l/m~i'. In the isotropic case, the odd functions 
S2n+i vanish, while for S2n dimensionality considerations give 

52„(r) = V"^'"-R2„(Ar,TOr), (6) 

where i?2n are some functions of dimensionless parameters. In principle, they can be calculated within the ordinary 
perturbation theory (i.e., as series in 50)1 but this is not useful for studying inertial-range behavior: the coefficients 
are singular in the limits Ar ^ 00 and/or mr 0, which compensate the smallness of 170 (assumed in ordinary 
perturbation theory), and in order to find correct IR behavior we have to sum the entire series. The desired summation 
can be accomplished using the field theoretic renormalization group (RG) and operator product expansion (OPE); 
see Refs. |2|-0. 

The RG analysis consists of two main stages. On the first stage, the multiplicative renormalizability of the model is 
demonstrated and the differential RG equations for its correlation functions are obtained. The asymptotic behavior 
of the functions like (||) for Ar > > 1 and any fixed mr is given by IR stable fixed points of the RG equations and has 
the form 

52„(r)- V"^'" (A0"^"^2„(™r), Ar»l, (7) 

with certain, as yet unknown, "scaling functions" R2n{'mr). In the theory of critical phenomena |^,^, the quantity 
A[S'2n] = — 2n + 7„ is termed the "critical dimension," and the exponent 7„, the difference between the critical 
dimension A[S'2n] and the "canonical dimension" — 2ri, is called the "anomalous dimension." In the case at hand, the 
latter has an extremely simple form: 7„ = ne. Whatever be the functions i?„(TOr), the representation (^ implies the 
existence of a scaling (scale invariance) in the IR region (Ar >> 1, mr fixed) with definite critical dimensions of all "IR 
relevant" parameters, A[5'2n] = —In + ne, = —1, A„i = 1 and fixed "irrelevant" parameters Vq and A. This means 
that the structure functions (|^) scale as S2n A'^['^^"l52n upon the substitution to A^^to, r A^t. In general, 
the exponent A[5'2n] is replaced by the critical dimension of the corresponding correlation function. This dimension is 
calculated as a series in e, so that the exponent e plays the part analogous to that played by the parameter e = A — d 
in the RG theory of critical phenomena, while i is an analog of the correlation length i = rc] see p2| , p3[ . 

On the second stage, the small mr behavior of the functions R2n{fnr) is studied within the general representation 
(^ using the OPE. It shows that, in the limit mr — s- 0, the functions R2n{mr) have the asymptotic forms 

R2nimr) = ^^^("^0 (mr)^[^l, (8) 

F 

where Cp are coefficients regular in mr. In general, the summation is implied over all possible renormalized composite 



operators F allowed by the symmetry (more precisely, see Sec. VII); A[F] being their critical dimensions. Without 
loss of generality, it can be assumed that the expansion is made in irreducible tensors; in the isotropic case only scalar 
operators have nonvanishing mean values and contribute to the right hand side of Eq. (^). 

The peculiarity of the models describing turbulence is the existence of the so-called "dangerous" composite operators 



with negative critical dimensions |L^ p^ , p4 25 1. Their contributions into the OPE give rise to a singular behavior of 
the scaling functions for mr —^ 0, and the leading term is given by the operator with minimal A[i^]. The leading 
contributions to 5*2™ are determined by scalar gradients F„ — [diOdiO)" and have the form 

524^)-^o~"^"^'"'^(™0^^ (9) 

where the critical dimensions A„ of the operators Fn are given by 

A„ = -2n{n - l)e/{d + 2) + 0{e^) ^ -2n(n - l)e/d + 0{l/d^). (10) 



The expression ^ agrees with the results obtained earlier in using the zero-mode techniques; the 0{e^) contri- 
bution to A„ is obtained in p^ . 



3 



In the theory of turbulence, the singular m dependence of correlation functions with the exponents nonlinear in n is 
referred to as anomalous scaling, and A„ themselves are termed the anomalous exponents; see, e.g., The above 

discussion shows that the anomalous exponents (in the sense of the turbulence theory) are not simply related to the 
anomalous or critical dimensions (in the sense of the theory of phase transitions) of the structure functions themselves; 
they are determined by the critical dimensions of certain composite operators entering into the corresponding OPE.Q 

In a number of papers, e.g. |l4 1 , the artificial stirring force in Eq. (|^) was replaced by the term (hv), where h 



is a constant vector that determines the distinguished direction and therefore introduces large-scale anisotropy. The 
anisotropy gives rise to nonvanishing odd functions S2n+i- The critical dimensions of all composite operators remain 
unchanged, but the irreducible tensor operators acquire nonzero mean values and their contributions appear on the 
right hand side of Eq. (^) ; see Q . This is easily understood in the language of the zero-mode approach: the noise 
/ and the term (hv) do not affect the differential operators in the equations satisfied by the equal-time correlations 
functions; the zero modes (homogeneous solutions) coincide in the two cases, but in the latter case the modes with 
nontrivial angular dependence should be taken into account. 

The direct calculation to the order 0(e) has shown that the leading exponent associated with a given rank contribu- 
tion to Eq. (|l|) decreases monotonically with the rank jlj]. Hence, the leading term of the inertial-range behavior of an 
even structure function is determined by the same exponent (|To|), while the exponents related to the tensor operators 
determine only subleading corrections. A similar hierarchy has been established recently in Ref. (see also [Q) 
for the magnetic field advected passively by the rapid-change velocity in the presence of a constant background field, 
and in within the context of the Navier-Stokes turbulence. 

Below we shall take the velocity statistics to be anisotropic also at small scales. We replace the ordinary transverse 
projector in Eq. (||) with the general transverse structure that possesses the uniaxial anisotropy: 

T,,(k) = a(V')Py (k) + 6(i^)n,(k)n,(k). (11) 

Here the unit vector n determines the distinguished direction (n^ = 1), ni(k) = P,;j(k)nj, and ip is the angle between 
the vectors k and n, so that (nk) = k cos ip [note that (fik) = 0] . The scalar functions can be decomposed the 
Gegenbauer polynomials (the d-dimensional generalization of the Legendre polynomials, see Ref. ]37|]): 

oo 



a(V') =^aiP2i(cosV'), 6(^) =^6,P2i(cosV') (12) 



1=0 1=0 

(we shall see later that odd polynomials do not affect the scaling behavior) . The positivity of the correlator (|^) leads 
to the conditions 

a(V') > 0, a(V') + 6(V') sin^ V > 0. (13) 

In practical calculations, we shall mostly confine ourselves with the special case 

T,,(k) = {l + pi cos^ ^)Py(k) + p2M^)hj{k). (14) 

Then the inequalities ( p^ reduce to pi,2 > —1. We shall see that this case represents nicely all the main features of 
the general model ( pi] ) 



We note that the quantities (|l l|) , (14) possess the symmetry n — n. The anisotropy makes it possible to introduce 
mixed correlator (v/) cx n6(t — t') C'{r/£) with some function C'{r/l) analogous to C{r/i) from Eq. (§). This violates 
the evenness in n and gives rise to nonvanishing o dd fu nctions 6*2,1+1. However, this leads to no serious alterations in 



the RG analysis; we shall discuss this case in Sec. VII, and for now we assume (v/) = 0. 

In a number of papers, e.g., |^0|-Q, the RG techniques were apphed to the anisotropically driven Navier-Stokes 
equation, including passive advection and magnetic turbulence, with the expression (|l4| ) entering into the stirring 
force correlator. The detailed account can be found in Ref. |2^, where some errors of the previous treatments are also 
corrected. However, these studies have up to now been limited to the first stage, i.e., investigation of the existence and 
stability of the fixed points and calculation of the critical dimensions of basic quantities. Calculation of the anomalous 
exponents in those models remains an open problem. 



^ The OPE and the concept of dangerous operators in the stochastic hydrodynamics were introduced and investigated in 
detail in |Q; see also |^,^. For the Kraichnan model, the relationship between the anomalous exponents and dimensions of 
composite operators was anticipated in within certain phenomenological formulation of the OPE, the so-called "additive 

fusion rules," typical to the models with multifractal behavior |^ . A similar picture has been discussed in | p8|j29| in connection 
with the Burgers turbulence and growth phenomena. 
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III. FIELD THEORETIC FORMULATION AND THE DYSON-WYLD EQUATIONS 

The stochastic problem (||)-(||) is equivalent to the field theoretic model of the set of three fields <i> = {9\ 6, v} with 
action functional 

S'($) = e'Dg9'/2 + 0' [^dt - (vO) + i^oA] 9 - (15) 

The first four terms in Eq. ( p^ represent the Martin-Siggia-Rose-type action for the stochastic problem (|^), (||) 
at fixed V (see, e.g., 12^,^), and the last term represents the Gaussian averaging over v. Here De and are the 
correlators (|^) and (^), respectively, the required integrations over x = (i,x) and summations over the vector indices 
are implied. 

The formulation ( |l5| ) means that statistical averages of random quantities in stochastic problem (|l|)-(||) can be 
represented as functional averages with the weight expS'(<i>), so that the generating functionals of total [G(j4)] and 
connected [VF(A)] Green functions of the problem are given by the functional integral 

G{A) = expM^(A) = j P$exp[S'($) + A$] (16) 

with arbitrary sources A = ,A^, A^ in the linear form 

A^= I dx[A^\x)9'{x)+A%x)9{x)+AJ{x)v,{x)]. 



The model ( |l^ ) corresponds to a standard Feynman diagrammatic technique with the triple vertex —9'{vd)9 = 
9'VjVj9 with vertex factor (in the momentum-frequency representation) 

VG- = -ifc„ (17) 

where k is the momentum flowing into the vertex via the field 9. The bare propagators in the momentum-frequency 
representation have the form 

{99')o = {9'9)* = {-tLo + iyokY\ {99)o ^ C{k) {lu^ + i^^k^', {9'9')o^0, (18) 

where C{k) is the Fourier transform of the function C(r/i) from Eg . (|^) and the bare propagator (vv)o = (vv) is 
given by Eq. (^ with the transverse projector from Eqs. ( pT| ) or ([l^). 

The pair correlation functions ($$) of the multicomponent field $ = {9',9,v} satisfy standard Dyson equation, 
which in the component notation reduces to the system of two equations, cf. [|l[ 

G'\cu, k) = -iuj + - Ee/e(cj, k), (19a) 

i?(w,k) = |G(c^,k)p[C(A:) + Ee'e'(cc',k)], (19b) 

where G'(tJ,k) = {99') and D{ijj,'k) = {99) are the exact response function and pair correlator, respectively, and S^'e, 
Tigigi are self-energy operators represented by the corresponding 1-irreducible diagrams; all the other functions 
in the model (|l^) vanish identically. 

The feature characteristic of the models like ( p^ ) is that all the skeleton multiloop diagrams entering into the 
self-energy operators E^'e, T^qiqi contain effectively closed circuits of retarded propagators {99') (it is crucial here that 
the propagator (vv)o in Eq. (H) is proportional to the 5 function in time) and therefore vanish. 

Therefore the self-energy operators in ( pj| ) are given by the one-loop approximation exactly and have the form 

E..(.,k)= ^ --Dohk,J'^J^^^^,^^^^Gi.',^'), (20a) 

duj' f dq Tij{q) 



I],,.,(.,k)= ^ =^ofc.fc.y (^(^^^Jy|^i?(c.',q'), (20b) 

where q' = k — q. The solid lines in the diagrams denote the exact propagators {99') and {99) , the ends with a slash 
correspond to the field 9', and the ends without a slash correspond to 9; the dashed lines denote the bare propagator 
(^; the vertices correspond to the factor (17). 
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The integrations over to' in the right-hand sides of Eqs. give the equal-time response function G(q) = 

(l/27r) / duj' G{uj' ,q) and the equal-time pair correlator D{q) = (l/27r) / dco' D{uj' ,q); note that both the self-energy 
operators are in fact independent of uj. The only contribution to G{q) comes from the bare propagator {69')q from 
Eq. (|l8|), which in the t representation is discontinuous at coincident times. Since the correlator which enters 
into the one- loop diagram for Sg'e, is symmetric in t and t\ the response function must be defined at t t' by half 



the sum of the limits. This is equivalent to the convention G(q) = (l/27r) / duj' {—iuj' 



= 1/2 and gives 



Ee'0(w,k) 



Dq kik. 



dq 



T^M) 



{2nY {q^ + m2)'i/2+e/2 



(21) 



The integration of Eq. ( 19b ) over the frequency w gives a closed equation for the equal-time correlator. Using Eq. 
(pl|) it can be written in the form 



2iyok^D(]i.) = C{k) 



Dq hkj 



dq 



^(q')-^(k) 



(22) 



Equation (|2^) can also be rewritten as a partial differential equation for the pair correlator in the coordinate 
representation, D{r) = {d{t,x.)6{t,x + r)) [we use the same notation D for the coordinate function and its Fourier 
transform]. Noting that the integral in Eq. (p2[) involves convolutions of the functions -D(k) and Dq T^ij (q) /(^^ + 



(23) 



2^d/2+e/2^ and replacing the momenta by the corresponding derivatives, iki —> di and so on, we obtain: 

2i^od^D{r) + C{r/£) + Dq S^j{r) d^dj D{r) = 0, 



where the "effective eddy diffusivity" is given by 

dq 



5^,(r) 



{2TTy (q2 + rn'^)d/2+e/2 



1 — exp (iq • r) 



(24) 



(note that diSij ~ 0). For < e < 2, equations (g^)-(^J) allow for the limit m ^ 0: the possible IR divergence of 
the integrals at q = is suppressed by the vanishing of the expressions in the square brackets. For the isotropic case 
(i.e., after the substitution Tij — > Pij) Eq. ( p3| ) coincides (up to the notation) with the well-known equation for the 
equal-time pair correlator in the model 

Equation ( pl| ) will be used in the next section for the exact calculation of the RG fimctions; solution of Eq. ( |25 
will be discussed in Sec. VIII in detail. 



IV. RENORMALIZATION, RG FUNCTIONS, AND RG EQUATIONS 

The analysis of the UV divergences is based on the analysis of canonical dimensions. Dynamical models of the type 
(p^), in contrast to static models, have two scales, so that the canonical dimension of some quantity F (a field or a 
parameter in the action functional) is described by two numbers, the momentum dimension dp and the frequency 
dimension d'^. They are determined so that [F] ~ [L]~'^i'[T]~'^p , where L is the length scale and T is the time scale. 
The dimensions are found from the obvious normalization conditions = —d^ = 1, d'^ — d"^ = 0, d^ = df = 0, 
= —df — 1, and from the requirement that each term of the action functional be dimensionless (with respect to the 
momentum and frequency dimensions separately). Then, based on dp and dp, one can introduce the total canonical 
dimension dp ^ dp + 2dp (in the free theory, 9t (x A), which plays in the theory of rcnormalization of dynamical 
models the same role as the conventional (momentum) dimension docs in static problems. 

The dimensions for the model ( [l5| ) are given in Table |, including the parameters which will be introduced later 
on. ^From Table || it follows that the model is logarithmic (the coupling constant go is dimensionless) at £ = 0, so 
that the UV divergences have the form of the poles in e in the Green functions. 

The total canonical dimension of an arbitrary 1-irreducible Green function F = ($ • • • is given by the relation 

dr=d^ + 2d^ =d + 2- N^d<j>, (25) 

where N,^, = {Nqi , Nq, Ny,} are the numbers of corresponding fields entering into the function F, and the summation 
over all types of the fields is implied. The total dimension dr is the formal index of the UV divergence. Superficial 
UV divergences, whose removal requires counterterms, can be present only in those functions F for which dr is a 
non-negative integer. 

Analysis of the divergences should be based on the following auxiliary considerations: 
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(i) From the explicit form of the vertex and bare propagators in the model ( |l5| ) it follows that Ngi — Ng = 2Nq 
for any 1-irreducible Green function, where Nq > is the total number of bare propagators {00) entering into the 
function (obviously, no diagrams with A'o < can be constructed). Therefore, the difference Ngi — Ng is an even 
non-negative integer for any nonvanishing function. 

(ii) If for some reason a number of external momenta occurs as an overall factor in all the diagrams of a given 
Green function, the real index of divergence d'^ is smaller than dr by the corresponding number (the Green function 



requires counterterms only if dp is a non-negative integer). In the model (15), the derivative d at the vertex 0'{vd)0 
can be moved onto the field 0' by virtue of the transversality of the field v. Therefore, in any 1-irreducible diagram 
it is always possible to move the derivative onto any of the external "tails" or 0' , which decreases the real index of 
divergence: d'^ = dr — Ng — Ngi . The fields 0, 0' enter into the counterterms only in the form of derivatives dO, dO'. 

From the dimensions in Table | we find dr = d + 2 - + Ng - {d + l)Ng, and d'j. = {d + 2){1 - Ng^) - N^. It 
then follows that for any d, superficial divergences can exist only in the 1-irreduciblc functions {0' . . . 0) i--^- with 
Ngi = 1 and arbitrary value of Ng, for which dr — 2, d'^ = 0. However, all the functions with Ng > Ngi vanish (see 
above) and obviously do not require counterterms. We are left with the only superficially divergent function {9'0)i-ii-; 
the corresponding counterterms must contain two symbols d and in the isotropic case reduce to the only structure 
0' AO. In the presence of anisotropy, it is necessary to also introduce new counterterm of the form 0'{nd)^0, which 
is absent in the unrenormalized action functional (|l5|). Therefore, the model ( |l5| ) in its original formulation is not 
multiplicativcly rcnormalizable, and in order to use the standard RG techniques it is necessary to extend the model 
by adding the new contribution to the unrenormalized action: 

S'($) = 0'DgO'/2 + 0' [-dt - (v9) + uqA + aoMrid)^] S - vD-^v/2. (26) 

Here ao is a new dimensionless unrenormalized parameter. The stability of the system implies the positivity of the 
total viscous contribution vok^ + Q;o^'o(nk)^, which leads to the inequality ao > —1. Its real ("physical") value is 
zero, but this fact does not hinder the use of the RG techniques, in which it is first assumed to be arbitrary and the 
equality ao = is imposed as the initial condition in solving the equations for invariant variables (see Sec. [V|). Below 
we shall see that the zero value of ao corresponds to certain nonzero value of its renormalized analog, which can be 
found explicitly. 

For the action (Eq), the nontrivial bare propagators in (nS) are replaced with 



{00% = {O'er, = (-i^ + '^ok' + aoMnkf)-' , {00)o = , . ^ ^^-^ . (27) 

\ - luj + i^ok"^ + aoi^o{ni^) \ 

After the extension, the model has become multiplicativcly rcnormalizable: inclusion of the counterterms is repro- 
duced by the inclusion of two independent renormalization constants Zi^2 as coefficients in front of the counterterms. 
This leads to the renormalized action of the form 

SR{<i>) = 0'DgO'/2 + 0' [-dt - (vd) + vZiA + avZ2{ndf] - vD-iv/2, (28) 

or, equivalently, to the multiplicative renormalization of the parameters vo, go and ag in the action functional (p6|): 

vo^vZy, go^giJL^Zg, ao^aZ^. (29) 

Here [i is the reference mass in the minimal subtraction (MS) scheme, which we always use in what follows, 
g and a are renormalized counterparts of the corresponding bare parameters, and the renormalization constants 
Z = Z((/,a; /9i,2, £, c') depend only on the dimensionless parameters. The correlator (^) in ( p8[ ) is expressed in 
renormalized variables using Eqs. (p9|). The comparison of Eqs. (p^), (p^), and ( |29| ) leads to the relations 

Z\ = Zi,, Zi = ZaZ^, Zg = Z^ ^. (30) 

The last relation in ( pO[ ) results from the absence of renormalization of the contribution with Do, so that Do = gova = 
g^ji^v', see (Q). No renormalization of the fields, anisotropy parameters and the "mass" m is required, i.e., Z$ — 1 for 
all and so on. 

The relation S{'^, gq) = Sb.{^, e, /i) (where cq = {vo, go, Q^o} is the complete set of bare parameters, and e = {j^, g, a} 
is the set of their renormalized counterparts) for the generating functional VK(A) in Eq. ( |l6|) yields W{A,eo) = 
Wfi{A, e, /i). We use to denote the differential operation fj,d^ at fixed cq and operate on both sides of this equation 
with it. This gives the basic RG differential equation: 

VRGWRiA,e,fi)=0. (31) 
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Here Vug is the operation P., expressed in the renormaHzed variables: 



(32) 



where we have written "Dx = xdx for any variable x, and the RG functions (the (3 functions and the anomalous 
dimension 7) are defined as 



for any renormalization constant Zp and 

I3g{g, a) EE Vf.g = g (-e - 7^) = g (-e + 7^) = ^ (_£ + ^^), 



(33a) 



(33b) 



I3a[g,a) = V^a = -aja = a(7i - 72)- 



(33c) 



The relation between (3g and j^, in Eq. ( p3b| ) results from the definitions and the last relation in (|30|). 

Now let us turn to the calculation of the renormalization constants Zi_2 in the MS scheme. They are determined 
by the requirement that the Green function G = (d'O), when expressed in renormalized variables, be UV finite, i.e., 
have no poles in e. It satisfies the Dyson equation of the form 



G ^(w,k) = -iw + z^ofc^ + aoJ^o(nk)^ - Ee/e(w,k), 



(34) 



which is obtained from Eq. (19a) with the obvious substitution i/qA;^ v^k^ + aofofnk)^; the self-energy operator 
Sg/g is given by the same expression (^ij). With Ty from ( p^ ) the integration in Eq. (|l]) yields 



Se'e(fc) 



Da J{m) 



2d{d + 2) 
where we have written 

J(m) 



(d - l){d + 2)fc2 + pi(d + l)fc2 + - 2pi(nk)2 + (^2 _ 2)p2(nk)^ 



dq 



r{e/2)m- 



(2^)^ (g2 + TO2)<i/2+^/2 (47r)''/2r(d/2 + e/2)' 
In deriving Eqs. (|35|), (^, we used the relations 



rfq/(9)^ 



d 



dqfiq), J dcifiq) 



Qi^jqilp _ S^jSip + SipSij + 5ii5p 



d{d + 2) 



dqfiq). 



(35) 



(36) 



(37) 



The renormalization constants Z1.2 are found from the requirement that the function {9' 9) expressed in renormalized 
variables be finite for e ^ 0. This requirement determines Zi^2 up to an UV finite contribution; the latter is fixed 
by the choice of a renormalization scheme. In the MS scheme all renormalization constants have the form "1 + only 
poles in e." We substitute Eqs. (p9|), (30) into Eqs. (|3^), ( ^5|) and choose Z12 to cancel the pole in e in the integral 
J{'m). This gives: 



gCd 



2d{d + 2)e 



{d-l){d+2)+ pi{d+l) + p2 



Z. 



1 



gCd 



2d{d + 2)ae 



-2pi + p2(rf'-2) 



(38) 



where we have written Cd = Sd/{2Tr'^) and Sd = 2t:'^/'^ lT{d/ 2) is the surface area of the unit sphere in d-dimensional 
space. 



For the anomalous dimension 71(17) = P^lnZi = j3gdg\nZi from the relations (33b) and ( |3q ) one obtains (note 
that Zi is independent of a) 



71(5) 



—sVg In Zi, 
l-Vg\nZ^ 



gCd 



2d{d + 2) 



(d- l)(d + 2)+pi(d+l) + p2 



(39) 



and for 72 (5, a) = T>^ In Z2 — iPgdg + Pada) In Z2 one has 



72(5,0) = 



[(-e + 7i)P<,+7i)I?, 



InZa 



-eP„ lnZ2 



gCd 



1 + 2?„ In Z2 



l + X'„lnZ2 2d{d + 2)a 



[-2pi+p2(d'-2)] 



(40) 
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[note that {Vg + 'Da)lTi Z2 = 0]. The cancellation of the poles in e in Eqs. (|39|), ( ^0| ) is a consequence of the UV 
finiteness of the anomalous dimensions jf; their independence of £ is a property of the MS scheme. Note also that 
the expressions (|38|)-(^0|) are exact, i.e., have no corrections of order and higher; this is a consequence of the fact 
that the one-loop approximation (|2l] ) for the self-energy operator is exact. 

The fixed points of the RG equations are determined from the requirement that all the beta functions of the model 
vanish. In our model the coordinates , a* of the fixed points are found from the equations 

Pg{9*,a*) = f3a{g*,a^) ^0, (41) 

with the beta functions from Eqs. (|3^). The type of the fixed point is determined by the eigenvalues of the matrix 
= {flik — dPi/dgk}, where (3i denotes the full set of the beta functions and g^ is the full set of charges. The IR 
asymptotic behavior is governed by the IR stable fixed points, i.e., those for which all the eigenvalues are positive. 



From the explicit expressions (39), (^0|) it then follows that the RG equations of the model have the only IR stable 
fixed point with the coordinates 

^ ^ 2did + 2)e -2pi+p2(rf'-2) 

^* " {d-l)(d + 2) + pi{d+l)+p2' * {d-l){d + 2) + pi{d+l) + p2- ^ ' 

For this point, both the eigenvalues of the matrix Q. equal to £; the values "ft — ~ it = ^ ^''''^ also found exactly 
from Eqs. (|3^ ) [here and below, 7^ = 'yrig*, «*)]• The fixed point ( ^2|) is degenerate in the sense that its coordinates 
depend continuously on the anisotropy parameters Pi,2-0 



V. SOLUTION OF THE RG EQUATIONS. INVARIANT VARIABLES 



The solution of the RG equations is discussed in Refs. |22|-|26|| in detail; below we confine ourselves to only the 
information we need. 

Consider the solution of the RG equation on the example of the even different-time structure functions 

52„(r,T) = ([0(i,x)-^?(t',x')]'"), r = x-x', t = t - t' . (43) 



It satisfies the RG equation 'DB.GS2n — with the operator Vrg from Eq. (^2|) (this fact does not follow automatically 
from Eq. ( ^l| ) and will be justified below in Sec. VI). 

In renormalized variables, dimensionality considerations give 

52„(r, t) = iy-"r^''R2n{tJ.r, Tv/r^,mr, g, a), (44) 

where i?2n is a function of completely dimensionless arguments (the dependence on d, e, pi^2 and the angle between 
the vectors r and n is also implied). ^From the RG equation the identical representation follows, 

52„(r, r) = (j?)-'V2"i?2n(l, TP/r\mr, g, a), (45) 

where the invariant variables e = e{pr, e) satisfy the equation Vjice = and the normalization conditions e = e at 
pr = I (we recall that e = {v, g,a,m} denotes the full set of renormalized parameters). The identity m = m is a 
consequence of the absence of Vm in the operator Vrg owing to the fact that m is not renormalized. Equation ( |45| ) 
is valid because both sides of it satisfy the RG equation and coincide for pr = 1 owing to the normalization of the 
invariant variables. The relation between the bare and invariant variables has the form 

iyo^i>Z^{g), go ^ gr"" Zg{g), ao = dZa{g,a). (46) 

Equation ( ^ ) determines implicitly the invariant variables as functions of the bare parameters; it is valid because both 
sides of it satisfy the RG equation, and because Eq. (^^ aX pr = 1 coincides with ( p9| ) owing to the normalization 
conditions. 



^ Formally, pi,2 can be treated as the additional coupling constants. The corresponding beta functions /3i,2 = Cjipi.a vanish 
identically owing to the fact that pi,2 are not renormalized. Therefore the equations /3i,2 = give no additional constraints on 
the values of the parameters g,a at the fixed point. 
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In general, the large /ir behavior of the invariant variables is governed by the IR stable fixed point: 5 — > 5*, a — > a, 
for /zr — > cxD. However, in multicharge problems one has to take into account that even when the IR point exists, 
not every phase trajectory (i.e., solution of Eq. (^)) reaches it in the limit ^ 00. It may first pass outside the 
natural region of stability [in our case, g > 0, a > —1] or go to infinity within this region. Fortunately, in our case 
the constants Zp entering into Eq. ( |4^ ) are known exactly from Eqs. (|3^), and it is readily checked that the RG flow 
indeed reaches the fixed point (|42| ) for any initial conditions go > 0, ao > —1, including the physical case ao = 0. 
Furthermore, the large fir behavior of the invariant viscosity P is also found explicitly from Eq. ( ^6| ) and the last 
relation in (|30|): D = D^r^ /g —> Dqt'^ /g^, (we recall that Dq — govo). Then for fxr 00 and any fixed rnr we obtain 

S2n{r,T) = D,^r^^^-'^g^^2n{rD„r^\mr), (47) 

where 

(.27i{DoTr^\mr) = R2n{l, DoTr^\mr, g^,a^), (48) 

and At = — 2 + 7* = — 2 + e is the critical dimension of time. The dependence of the scaling function ^2ri on its 
arguments is not determined by the RG equation ( pl| ) itself. For the equal-time structure function (|^), the first 
argument of ^2n in the representation ( ^8|) is absent: 

52„(r) = i?„-"r"(2-e)gn ^2„(mr), (49) 

where the definition of ^2n is obvious from (^). It is noteworthy that Eqs. (^^-(^9|) prove the independence of the 
structure functions in the IR range (large jir and any mr) of the viscosity coefficient or, equivalently of the UV scale: 
the parameters go and 1^0 enter into Eq. ( p7| ) only in the form of the combination Dq — goi/Q. A similar property was 
established in Ref. |35j for the stirred Navier-Stokes equation and is related to the Second Kolmogorov hypothesis; 
see also the discussion in ||25| , ^ . 

Now let us turn to the general case. Let F{r, r) be some multiplicatively renormalized quantity (say, a correlation 
function involving composite operators), i.e., F = ZpFn with certain renormalization constant Zp. It satisfies the 
RG equation of the form [I?_rg + Jf]Fii — with jp from ( ^3a|) . Dimensionality considerations give 

Fnir^T) = v'^^^r-'^''Rp{fjLr,Tv/r^,mr,g,a), (50) 



where and dp are the frequency and total canonical dimensions of F (see Sec. IV) and Rp is a function of 
dimensionless arguments. The analog of Eq. ( ^ ) has the form 

F(r,T) = Zp{g,a)FR = Zp{g,a) {D^^^r-''^ Rp{l,TD/r^ ,mr,g,a). (51) 

In the large fir limit, one has Zp{g,a) ~ const (Ar) ; see, e.g., [P6| . The UV scale appears in this relation from 
Eq. (^). Then in the IR range (Ar ~ fir large, mr arbitrary) Eq. (^1[) takes on the form 

F(r,r) ~ const A-''^-Do^r-^['^leF(-DoTr^Smr). (52) 

Here 

A[F]=Ap = d''p-Atd^p + j*p, At^-2 + e (53) 

is the critical dimension of the function F and the scaling function ^p is related to Rp as in Eq. (|47|). For nontrivial 
7^, the function F in the IR range retains the dependence on A or, equivalently, on vq. 



VI. RENORMALIZATION AND CRITICAL DIMENSIONS OF COMPOSITE OPERATORS 

Any local (unless stated to be otherwise) monomial or polynomial constructed of primary fields and their deriva- 
tives at a single spacetime point x = (t,x) is termed a composite operator. Examples are 6^{x), [di9{x)did{x)]^ , 
di9{x)dj9{x), e'{x)Vt9{x) and so on. 

Since the arguments of the fields coincide, correlation functions with such operators contain additional UV diver- 
gences, which are removed by additional renormalization procedure; see, e.g., [ p2| , p3[ . For the renormalized correlation 
functions standard RG equations are obtained, which describe IR scaling of certain "basis" operators F with definite 
critical dimensions Ap = A[F]. Due to the renormalization, A[F\ does not coincide in general with the naive sum of 
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critical dimensions of the fields and derivatives entering into F. As a rule, composite operators mix in renormaliza- 
tion, i.e., an UV finite renormalized operator has the form F^ — F+ counterterms, where the contribution of the 
counterterms is a linear combination of F itself and, possibly, other unrenormalized operators which "admix" to F in 
renormalization. 

Let F = {Fa} be a closed set, all of whose monomials mix only with each other in renormalization. The renormal- 
ization matrix Zp = {Zap} and the matrix of anomalous dimensions ^p = {70/3} for this set are given by 

Fa=Yl ZapF^', IF = Zp^V^Zp, (54) 

13 

and the corresponding matrix of critical dimensions Ai^ = {Aq/j} is given by Eq. (|53|), in which and df^ are 
understood as the diagonal matrices of canonical dimensions of the operators in question (with the diagonal elements 
equal to sums of corresponding dimensions of all fields and derivatives constituting F) and 7^ = 7_F(g*,a*) is the 
matrix (|5^) at the fixed point (|2|) . 

Critical dimensions of the set F = {Fa} are given by the eigenvalues of the matrix Ap. The "basis" operators that 
possess definite critical dimensions have the form 

fj 

where the matrix Up — {Ua/3} is such that A^ = UpApUp^ is diagonal. 

In general, counterterms to a given operator F are determined by all possible 1-irreducible Green functions with 
one operator F and arbitrary number of primary fields, T — {F{x)^{xi) . . . ^{x2))i-h-. The total canonical dimension 
(formal index of UV divergence) for such function is given by 

dr ~ dp — N^d,f,, (56) 

with the summation over all types of fields entering into the function. For superficially divergent diagrams, dr is a 
non-negative integer, cf. Sec. |^ 

Let us consider operators of the form 9^{x) with the canonical dimension dp = —N, entering into the structure 
functions (|^) . From Table | in Sec. IV and Eq. ( |5^ we obtain dr = —N + Ng — — {d+ l)Ng> , and from the analysis 



of the diagrams it follows that the total number of the fields 9 entering into the function F can never exceed the 
number of the fields 9 in the operator 6^ itself, i.e., Ng < N (cf. item (i) in Sec. 0). Therefore, the divergence can 
only exist in the functions with = Ng> = 0, and arbitrary value of = Ng, for which the formal index vanishes, 
dr — 0. However, at least one of Ng external "tails" of the field 9 is attached to a vertex 6'{vd)6 (it is impossible 
to construct nontrivial, superficially divergent diagram of the desired type with all the external tails attached to the 
vertex F), at least one derivative d appears as an extra factor in the diagram, and, consequently, the real index of 
divergence dp is necessarily negative. 

This means that the operator 6^ requires no counterterms at all, i.e., it is in fact UV finite, 0^ = Z [9^]^ with 
Z ^ 1. It then follows that the critical dimension of 9^ (x) is simply given by the expression (|53| ) with no correction 
from "fp and is therefore reduced to the sum of the critical dimensions of the factors: 

A[9^]^NA[9]=N{~l + e/2). (57) 

Since the structure functions (^) or ( ^ ) are linear combinations of pair correlators involving the operators 9-^ , equation 
(57) shows that they indeed satisfy the RG equation of the form ([3TI), discussed in Sec. 0. We stress that the relation 
(57) was not clear a priori; in particular, it is violated if the velocity field becomes nonsolenoidal [ p^ . 

In the following, an important role will be also played by the tensor composite operators di-^9 ■ ■ ■ di^9 {di9di9)^ 
constructed solely of the scalar gradients. It is convenient to deal with the scalar operators obtained by contracting 
the tensors with the appropriate number of the vectors n, 

F[N,p] = [{n^)9]P{^^9^^9)'', N=2n+p. (58) 

Their canonical dimensions depend only on the total number of the fields 9 and have the form dp — 0, dp — ~N. 

In this case, from Table | and Eq. (|5^ ) we obtain dr = Ng — N^ — {d + l)Ngi, with the necessary condition 
Ng < N, which follows from the structure of the diagrams. It is also clear from the analysis of the diagrams that the 
counterterms to these operators can involve the fields 9, 9' only in the form of derivatives, 89, 89' , so that the real 
index of divergence has the form d'r = dr — Ng — Ng' — —N^ — {d + 2)Ngi . It then follows that superficial divergences 
can exist only in the Green functions with N^ = Ngi = and any Ng < N, and that the corresponding operator 
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counterterms reduce to the form F[N' ,p'] with A^' < N. Therefore, the operators ( |5q ) can mix only with each other 
in renormahzation, and the corresponding infinite renormalization matrix 

F[N,p] = ^[N,p] [N',p'] F^'W^p'] (59) 

is in fact block-triangular, i.e., .^[Ar.p] [N',p'] — for N' > N . It is then obvious that the critical dimensions associated 
with the operators F[N^p\ are completely determined by the eigenvalues of the finite subblocks with N' — N . In 
the following, we shall not be interested in the precise form of the basis operators (|55|), we rather shall be interested 
in the anomalous dimensions themselves. Therefore, we can neglect all the elements of the matrix ( ^9|) other than 

Z[N,p] [N,p']- 

In the isotropic case, as well as in the presence of large-scale anisotropy, the elements ■^[Ar,p] [N,p'] vanish for p < p' , 
and the block ^[w.p] [N,p'] is triangular along with the corresponding blocks of the matrices Up and Ap from Eqs. ([5^), 
(^). In the isotropic case it can be diagonahzed by changing to irreducible operators (scalars, vectors, and traceless 
tensors), but even for nonzero imposed gradient its eigenvalues are the same as in the isotropic case. Therefore, 
the inclusion of large-scale anisotropy does not affect critical dimensions of the operators (|5|) ; see . In the case 
of small-scale anisotropy, the operators with different values of p mix heavily in renormalization, and the matrix 
Z[N,p] [N.p'] is neither diagonal nor triangular here. 

Now let us turn to the calculation of the renormalization constants ■Z^[Ar,p] [N.p'] in ths one- loop approximation. Let 
T(x; 0) be the generating functional of the 1-irreducible Green functions with one composite operator F[N,p] from 
Eq. ( p8|) and any number of fields 9. Here x = {t, x) is the argument of the operator and 6 is the functional argument, 
the "classical counterpart" of the random field 9. We are interested in the A^-th term of the expansion of r(x; 9) in 
9, which we denote Tn(x;9); it has the form 

^n{x;9) = l^Jdxi--- J dxN9{xi)---9ixN){F[N,p]{x)9{xi)---9ixN))i-ir- (60) 
In the one-loop approximation the function (^) is represented diagrammatically in the following manner: 

rN^F[N,p] + ^ ■ (61) 

Here the solid lines denote the bare propagator {99') o from Eq. (p7|), the ends with a slash correspond to the field 9' , 
and the ends without a slash correspond to 9; the dashed line denotes the bare propagator (^; the vertices correspond 
to the factor (|l^. 

The first term is the "tree" approximation, and the black circle with two attached lines in the diagram denotes the 
variational derivative V{x; xi,X2) = S''F[N,p]/S9{xi)S9{x2)- It is convenient to represent it in the form 

V{x; XI, X2) = d^S{x - xi) dj5{x - X2) [(na)f (a^)"] , (62) 

where is a constant vector, which after the differentiation is substituted with di9{x). The diagram in Eq. ( |6l| ) is 
written analytically in the form 

dxi--- / dxiV{x; Xi, X2){9{xi)9' {x3))o{9{x2)9' {xi))o{vk{x3)vi{x^))odk9{x3)di9{x4), (63) 



with the bare propagators from Eqs. (||), (|27|); the derivatives appear from the ordinary vertex factors (p^). 

In order to find the renormalization constants, we need not the entire exact expression (|6^), rather we need its UV 
divergent part. The latter is proportional to a polynomial built of N factors 89 at a single spacetime point x. The 
needed N gradients has already been factored out from the expression (|3|): {N — 2) factors from the vertex (|6^ ) and 
2 factors from the ordinary vertices (|l^). Therefore, we can neglect the spacetime inhomogeneity of the gradients and 
replace them with the constant vectors a^. Expression (|6^) then can be written, up to an UV finite part, in the form 

d' 



daidaj 
where we have denoted 



(na)f(a^)" X^.^u, (64) 



Xtj,ki= dx3 dx4 di{9{x)9'{x3))odj{9{x)9'{x4)}o{vk{x3)vi{x4))o, (65) 
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or, in the momentum-frequency representation, after the integration over the frequency. 



X., 



Do 



Tfei(q) 



21^0 J (27r)'i (q2 + a(qn)2) (^k^ + m^)-d/2-e/2 
with Dq from Eq. (0) and Tki from Eq. (p^). The tensor Xij^ki can be decomposed in certain basic structures: 



(66) 



X 



ij, kl 



ij, kl ' 



(67) 



where are scalar coefficients and S"*^"^ are all possible fourth rank tensors constructed of the vector n and Kronecker 
delta symbols, and symmetric with respect to the permutations within the pairs {ij} and {kl}: 



^ij!ki = (n-iUkSji + runiSjk + njUiSik + njnkSii)/A, s'^^^'^.i = runjUkni. 
It is now convenient to change from the tensors (|66|), (|67|) to the scalar integrals defined by the relations 



:(2) 



(3) 



;'(4) 



A, 



ijkl 



l,---,6 



(68) 



(69) 



[it is obvious from Eqs. (|67|), ( |68| ) that the coefficients A2 and vanish identically by virtue of the relation 
9/21fc(q) = 0]. The UV divergence and the angular integration in Eqs. (|6^), (|6^ ) can be disentangled using the 
relations 



(70) 



where the integral J(m) from (36) contains a pole in e, and Hg are UV finite and e independent quantities given by 
only the angular integrations: 



Hi 



Hi 



1 + a cos^ 

Tu cos^ -0 
1 + a cos^ 



H. 



He, 



1 + a cos^ tp 



niHjTij cos^ ^0 
1 + a cos^ ip 



(71) 



where -i/i is the angle between the vectors q and n, so that (nq) ~ qcosip, and the quantities Tu and niUjTij with Ty 
from (^) depend only on ^ (the dependence on d and the anisotropy parameters pi.2 is of course also implied). The 
brackets {■ ■ ■)s denote averaging over the unit d-dimensional sphere, which for a function dependent only on ^ takes 
on the form 



{fms = 



Sd- 



Sd 



di, sin'^-2^/(0). 



with Sd = /T{d/2). Finally, the change of variables y = cos^ i]j in the integrals ( fzi] ) gives 



dy — — ; — = B{a, h) F{c, a;a + b; —a), 







{l + ayY 

where B is the Euler beta function and F is the hypergeometric function, see, e.g. 



(72) 



(73) 



This gives: 



Hi 
Ha 
Hi 



{d^l + P2)Fo+ (rf-l)pi-p2 



El 
d ' 



Fi 3Fo 

(1 + P2) Fo + (-1 + PI - 2p2) - + ip2- P.)-^y 



Fi r 

{d-l + p2) — + (d-l)pi-p2 
d I 



Fi 



d{d + 2 
3F2 



Fa = (l + P.)^ + (-l + Pi-2p2)^^(^ + (p2-pi) 



I5F3 



d(d + 2)(d + 4)' 



(74) 



13 



where we have denoted F„ = 1/2 + n; d/2 + n; —a). 

In principle, all the hypergeometric functions entering into Eqs. (|7^ ) can be reduced to the only function, for 
example F3, but the coefficients then become too complex. 

After quite lengthy but straightforward calculations, the quantity in Eq. (^) can be expressed through the coeffi- 
cients Bs, then through A^, and finally through Hs'. 



^^J{m)[Q^F[N,p^2] + Q2F [N, p] + F[N,p + 2] + Q^F[N ,p + A]^, 
where we have substituted back with the gradients di9{x), used the notation (pq), and denoted 



(75) 



Qi 



1) V 



p{p - 1) 
(d- 

p{p - f ) 
(d-l) 



H4 — 



dHf. — Ua 



(N-p) 



(d-l)(d+l) 



-JiN-p 



d-l){Hi 

Hi [-{N - 2) + p{2d + 3)] + 3^6 [{N - 2) - p{d + 2)] }, 



Q: 



(N-p) 



-{-{2N^2p + d- 3)Hi + H3 [{N -p){d + 3) + {d + 2){d - 3)] - 

p(d + 4)]}, 



{d~l){d 

+ Hi [{N - 2)(d + 3) - p{3d + 8)] + 2HG{d + 2) [-3(7V - 2 



{N~p)iN~p-2) 



(d- l)(d+l) 



H, 



{d + 2)iHs + Hi) + {d + 2){d + A) He 



(76) 



The renormalization constants -^[Af^p] iN,p'] can be found from the requirement that the function ( |6l| ) after the 
replacement F[N,p] —>■ F^'[N,p] with F^\N,p] from ( |59| ) be UV finite, i.e., have no poles in e, when expressed in 



renormalized variables using the formulas (pQl). In the MS scheme, the diagonal element has the form -^[Ar^p] [N,p'' 
only poles in e, while the nondiagonal ones contain only poles. Then from Eqs. (|6^), ( |70| ) and (^5|) one obtains 



= 1- 



gCd 
As 

gCd 

4e 

with coefficients Qi from Eq. (|7^ ) and Crf from 
^[w.p] [Ar",p"]^M^[iV",p"] [W',?)'] ^^ave the form 



Z[N,p][N,p] — 1 



^[Af,p][Af,p+2] 



Q2 



[N,p][N,p-2] 



^[N ,p][N ,p+i] 



4e 

gCd 

4e 



(77) 



The corresponding anomalous dimensions 1[n,p][N' .p'\ = 



i[N,p][N,p] — -gCdQ2/^, i[N,p][N,p-2] — -gCdQi/A, 
7[7v,p][Ar,p+2] = -gCdQa/A, 7[Af,p][JV,p+4] = -gCdQi/A. 



(78) 



In contrast to exact expressions (^8|), the quantities in Eqs. ([ttI), ( [ts]) have nontrivial corrections of order and 
higher. The matrix of critical dimensions (^3|) is given by 



A[Ar,p][Ar,p'] - 7V£/2 l[N,p][N,p']' 



(79) 



where the asterisk impHes the substitution ([4 21). The set of equations (|42|), (|74D , ([76|), (|78|), ( |79D gives the desired 
expression for the matrix of critical dimensions of the composite operators (p8|) in the one-loop approximation, i.e., 
in the first order in e (we recall that 5* = 0{e) and a* = 0(1)). In contrast to simpler expressions like (|5^), the 
dimensions in Eq. ( [79|) depend on the anisotropy parameters pi^2 and have nontrivial corrections of order and 
higher. 

In the one-loop approximation, the expression (75) contains only four terms, and therefore all the matrix elements 
7[Af.p][Af,p'] other than ( fz^ ) vanish. However, they become nontrivial in higher orders, so that the basis operators ([5^), 
associated with the eigenvalues of the matrix (|79|), are given by mixtures of the monomials (58) with all possible 
values of the index p, or, in other words, mixtures of the tensors belonging to different representations of the rotation 
group. 



14 



As already said above, the critical dimensions themselves are given by the eigenvalues of the matrix (]79|). One can 
check that for the isotropic case (pi,2 = 0), its elements with p' > p vanish, the matrix becomes triangular, and its 
eigenvalues are simply given by the diagonal elements A[A^,p] = A[jv,p][Ar,p]- They are found explicitly and have the 
form 

It is easily seen from Eq. (|80| ) that for fixed N and any d> 2, the dimension A[A^, p] decreases monotonically with p 
and reaches its minimum for the minimal possible value of p = pN, i.e., p^ = if A'^ is even and pn = 1 if iV is odd: 

A[N,p]> A[N,p'] if p>p'. (81a) 

Furthermore, this minimal value A[A^, pjv] decreases monotonically as N increases for odd and even values of N 
separately, i.e., 

> A[2n,0] > A[2n + 2,0], A[2n + 1, 1] > A[2n + 3, 1]. (81b) 

A similar hierarchy is demonstrated by the critical dimensions of certain tensor operators in the stirred Navier-Stokes 
turbulence; see Ref. and Sec. 2.3 of However, no clear hierarchy is demonstrated by neighboring even and 
odd dimensions: from the relations 

A[2n +1,1]- A[2n, 0] = "^^/^'^''^ A[2n + 2,0]- A[2n +1,1] = (81c) 

it follows that the inequahty A[2n + 1, 1] > A[2n + 2, 0] holds for any d> 2, while the relation A[2n, 0] > A[2n + 1, 1] 
holds only if n is sufficiently large, n > (d + 2)/4. | 

In what follows, we shall use the notation A[Af,p] for the eigenvalue of the matrix (|7^) which coincide with (|8^) for 
pi_2 = 0. Since the eigenvalues depend continuously on pi,2, this notation is unambiguous at least for small values of 

Pl,2- 

The dimension A[2,0] vanishes identically for any pi 2 and to all orders in e. Like in the isotropic model, this can 
be demonstrated using the Schwingcr equation of the form 

j [eix) exp Sr{^) + A$] /69'{x) = 0, (82) 

(in the general sense of the word, Schwinger equations are any relations stating that any functional integral of a total 



variational derivative is equal to zero; see, e.g., |22 2^). In (p2|), is the renormalized action (P8[), and the notation 
introduced in ( [T^ ) is used. Equation ( p^ ) can be rewritten in the form 

e'Dg9 - Vt[d^/2] + iyZiA[e^/2] + aiyZ2ind)^[0^/2] - iyZiF[2, 0] - aiyZ2F[2, 2]\ = -~Ae'SWRiA)/6Ag. (83) 



Here Dg is the correlator (||), (• • ■)a denotes the averaging with the weight exp[S'fl;(<l') + A^], Wr is determined by 
Eq. ( p^ ) with the replacement S — > Sr, and the argument x common to all the quantities in (^) is omitted. 

The quantity {F) a is the generating functional of the correlation functions with one operator F and any number 
of the primary fields <&, therefore the UV finiteness of the operator F is equivalent to the finiteness of the functional 
{F)a- The quantity in the right hand side of Eq. ( |83| ) is UV finite (a derivative of the renormalized functional with 
respect to finite argument), and so is the operator in the left hand side. Our operators F[2, 0], F[2, 2] do not admix in 
renormalization to 9'Dg9 (no needed diagrams can be constructed), and to the operators Vt[6'^/2] and A[6'^/2] (they 
have the form of total derivatives, and F[N,p] do not reduce to this form). On the other hand, all the operators in 
( ^3| ) other than F[N,p] do not admix to F[N,p], because the counterterms of the operators (|5|) can involve only 
operators of the same type; see above. Therefore, the operators F[N,p] entering into Eq. (^3|) are independent of the 



^The situation is different in the presence of the linear mean gradient: the first term Ne/2 in Eq. (pOl ) is then absent owing 
to the difference in canonical dimensions, and the complete hierarchy relations hold, A[2n, 0] > A[2n + 1, 1] > A[2n + 2, 0]; see 



0- 
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others, and so they must be UV finite separately: vZiF[2^ 0] + avZ2F[2, 2] = UV finite. Since the operator in (83) is 
UV finite, it coincides with its finite part, 

vZiF[2, 0] + avZ2F[2, 2] = vF^[2, 0] + auF'^[2, 2], 

which along with the relation (|5^) gives 

ZlZ[2,Q][2,Q] + OlZ2Z\^2,2][2.0] = li ^'i Z[2^o] [2,2] + "■^2-Z'[2,2] [2,2] = Oi, 

and therefore for the anomalous dimensions in the MS scheme one obtains 

7l + 7[2,0][2,0] + a7[2,2][2,0] = 0, 7[2,0][2,2] + "72 + a7[2,2][2,2] = 0. 



These relations can indeed be checked from the explicit first-order expressions ([r7|), (|78|). Bearing in mind that 
7* = 72 = £ (see Sec. IV), we conclude that among the four elements of the matrix 7^, only two, which we take to be 
^'[2 2] [2 0] ^'^'^ "^[2 2] [2 2]' independent. Then the matrix of critical dimensions (|8^) takes on the form 



M2,p][2,p'] 



Q^*7[2,2][2,o] 

'>'[2,2][2,0] 



a*e '^*7[2,2][2,2] 

'T[2,2][2,2] 



(84) 



It is then easily checked that the eigenvalue of the matrix (|84D, which is identified with A [2,0], does not involve 
unknown anomalous dimensions and vanishes identically, A [2,0] = 0, while the second one is represented as 

A[2, 2] = £ - a*7*2,2][2,0] + 7[*2,2][2,2]- 

Using the exphcit 0{e) expressions ([76|), (|78|), one obtains to the order 0[e): 

A[2, 2]/£ = 2 + |-(d - 2)d{d + 2){d + 4)Fo* - (d + 2)(d + 4)(2 + (d - 2)pi + dp2)Fl + 
+3(d + 4) (d - 2pi + 2dp2)F; + 15d(pi ~ p2)F*y [{d-l){d + A)[{d-l){d + 2) + {d+l)pi+ p2]Y (85) 

where F* = F(l, 1/2 + n; d/2 + n; —a*) with a* from Eq. ( |4^ ) and the hypergeometric function from Eq. ([73|). 

In Fig. 1, we present the levels of the dimension (|85| ) on the pi-p2 plane for d = 3. We note that the dependence 
on pi_2 is quite smooth, and that A[2, 2] remains positive on the whole of the pi^2 plane, i.e., the first of the hierarchy 
relations (|8|) remains valid also in the presence of anisotropy. A similar behavior takes place also for d — 2. 

For N > 2, the eigenvalues can be found analytically only within the expansion in pi_2- Below we give such 
expansions up to the terms 0(pi 2) for all the eigenvalues with iV < 4 (with the notation p^ = 2pi + dp2)- 



7*[2,2]/£ = 



4(d-2)(d+l)p3 



2(d+ 1) 



(d-l)(d + 2) (d-l)2(d + 2)2(d + 4) (d-l)3(d + 2)3(d + 4)(d + 6) 



-6d(d2 + 6d - 15)pi + (d^ + 6d* - 41d^ + A2d'^ + 16d + 24)pip2 + (-lid'' + 23d-'^ + 2d^ - 26d + 24)p^ 



(86a) 



7*[3,l]/£ = - 



(d + 4) 2(d-2)(d+l)p3 



(d + 2) 



(d-l)(d + 2)3 
2(d5 + 15d'' + 



2(d- 2)(d+ 1) 
(d- l)2(d + 2)e(d + 4)2 

^ f 204d^ + 16d - m)pl 



-d{-d^ + d^ + 80d'* + 296d3 + 64d^ - 48d + M)pip2 + (3d^ + 29d^ 



74d'' - 32d3 - 72d2 + 48d + 64)/?^ 



(86b) 



7*[3,3]/£: 



12(d-2)(d+l)2p3 



6(d- 2)(d+ 1) 



(d-l)(d+2) (d- l)2(d+2)3(d + 4) (d- l)3(d + 2)6(d + 4)2(d + 6) 
X 2(3d*' + 54d^ + 321d'' + 722d^ + 420d2 + 184d + %6)p\ 
+d{-d' - 16d^ - 49d5 + 158d'^ + 844d3 + 904d2 + 1120d + 640)pip2 
+ (13d^ + llld^ + 258d^ + 26d'' - 124d^ + ^Sd^ + 336d + 192)p2 



(86c) 
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^ ^ IdT^ - (d-l)Mrf + 2)^(^ + 4) ' ^^''^^ 

_ -d^-5d + 8 4(d-2)(d+l)(d^+3d-16)p3 2{d+l) 

^ ^ ' {d^ l){d + 2) {d- l)2(d + 2)2(d + 4)2 (d - l)3d(d + 2)3(d + 4)5(d + 6) ^ 

-2(2d^ + 47d^ + 422d^ + 1253d^ - 3684^^ - 23344d'' - 4032^^ + 57856^^ + 86016d - 73728)p? + 

+d{d^ - 28d^ - 75M^ - 4654d^ - 2712d^ + 40232d* + 42592^^ - 76928^^ - 202240d + 122880)pip2 - 

-d{8d^ + U9d^ + 699d^ - 1074^'^ - 10994rf^ - 2240<i^ + 29248d^ + 29696^^ - 49664d + 24576)p^ , (86e) 

12 , 24(rf-2)(d+l)p3 12(d-2) 

7 [4, 4J/£ = 



(d-l)(d + 2) (d- l)2(d + 2)(d + 4)2 (d- l)3(d + 2)4(d + 4)5(d + 6) 
X 2d(3d^ + 81d® + 823d^ + 4023^^^ + 10230rf^ + 14536rf^ + 11040d + 3456)pi 
+ {-d^° ~ 23d^ - 163d^ - 225d^ + 2040d^ + 9948d^ + 21240d'' + 25424d^ + 18112^^ + 8960rf + 3072)pip2 
+ {15d^ + 232d^ + 1309d^ + 3306d^ + 3782^^ + 780d* - 768^^ + 3712^^ + 6656<i + 3072)p^ 



(86f) 



These expressions illustrate two facts which seem to hold for all N: 

(i) The leading anisotropy correction is of order O(/0i,2) for p ^ and 0{pi 2) for p = 0, so that the dimensions 
7* [N, 0] are anisotropy independent in the linear approximation, and 

(ii) This leading contribution depends on pi 2 only through the combination p3 = 2pi + dp2- 
This conjecture is confirmed by the following expressions for iV = 6, 8 and p = 0: 

Vrfi nl/.-^^i^±^ 12(rf-2)2(d+l)(d2 + 14rf + 48)pi 

^ ~ ~(dT2) (d-l)2d(d + 2)4(d + 4)2 ' ^^^^^ 

nsnl/. "4(^ + 8) 24(d-2)2(d+l)(d2 + i8d + 80)p2 
^ IdT^^ (d-l)2d(d + 2)4(d + 4)2 ' 

The eigenvalues beyond the small pi_2 expansion have been obtained numerically. Some of them are presented in 
Figs. 2-5, namely, the dimensions A[n,p] for n = 3,4,5,6 vs pi for p2 — 0, vs pi = p2, and vs p2 for pi = 0. The 
main conclusion that can be drawn from these diagrams is that the hierarchy (p^) demonstrated by the dimensions 
for the isotropic case (pi_2 = 0) holds valid for all the values of the anisotropy parameters. 

VII. OPERATOR PRODUCT EXPANSION AND ANOMALOUS SCALING 

Representations (^9|) for any scalingfunctions ^{mr) describes the behavior of the correlation functions for Ar >> 1 
and any fixed value of mr; see Sec. |V|. The inertial range corresponds to the additional condition mr « 1. The 
form of the functions ^{mr) is not determined by the RG equations themselves; in analogy with the theory of critical 
phenomena, its behavior for mr — > is studied using the well-known Wilson operator product expansion (OPE); see, 
e.g., 

According to the OPE, the equal-time product Fi {x)F2 {x') of two renormalized operators for x = (x-|-x') /2 = const 
and r = x — x' ^ has the representation 



F,{x)F2(x')^Y.^F{v)F{t,^), 



where the functions Cp are coefficients regular in and F are all possible renormalized local composite operators 
allowed by symmetry (more precisely, see below). Without loss of generality, it can be assumed that the expansion 
is made in basis operators of the type (pq), i.e., those having definite critical dimensions Ap. The renormalized 
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correlator {Fi{x)F2{x')) is obtained by averaging Eq. dSSj ) with the weight expSa with the renormaUzed action (|2? 
the quantities (F) cx m^^ appear on the right hand side. 

From the operator product expansion ( p8| ) we therefore find the fohowing expression for the scahng function ^ (mr) 
in the representation ( ^9|) for the correlator {Fi{x)F2{x')): 

^{mr) ^^Af (mr)^^ , mr « I, (89) 
F 

with the coefficients Ap regular in (mr)'^. 

Now let us turn to the equal-time structure functions Sn from (H). From now on, we assume that the mixed 
correlator (v/) differs from zero (see Sec. ||); this does not affect the critical dimensions, but gives rise to nonvanishing 
odd structure functions. In general, the operators entering into the OPE are those which appear in the corresponding 
Taylor expansions, and also all possible operators that admix to them in renormalization p2| , p3| . The leading term of 
the Taylor expansion for the function Sn is obviously given by the operator F[N, N] from Eq. (^q); the renormalization 
gives rise to all the operators F[N',p\ with N' < N and all possible values of p. The operators with N' > N (whose 
contributions would be more important) do not appear in Eq. (|89|), because they do not enter into the Taylor 



expansion for Sn and do not admix in renormalization to the terms of the Taylor expansion; see Sec. VI. Therefore, 
combining the RG representation ( ^7|) with the OPE representation (^9|) gives the desired asymptotic expression for 
the structure function in the inertial range: 

^^(r) = ^/V^'d-/^) ^ J2{Cn',, (mr)^[~'-^l + ..•}. (90) 

N'<N p 

The second summation runs over all values of p, allowed for a given N'; C'n'.p arc numerical coefficients dependent on 
e, d, pi, 2 and the angle d between r and n. The dots stand for the contributions of the operators other than F[N,p], 
for example, d^9d^6] they give rise to the terms of order [mrY^'-"^'^^ and higher and will be neglected in what follows. 
Some remarks are now in order. 

(i) If the mixed correlator (v/) is absent, the odd structure functions vanish, while the contributions to even 
functions are given only by the operators with even values of N' . In the isotropic case (pi,2 = 0) only the contributions 
with p = Q survive; see In the presence of the anisotropy, pi_2 7^ Oi the operators with p 7^ acquire nonzero 
mean values, and their dimensions A[iV',p] also appear on the right hand side of Eq. (|90|). 

(ii) The leading term of the small mr behavior is obvi ously given by the contribution with the minimal possible 



value of A[7V',p]. Now we recall the hierarchy relations (31a), ( Sib ), which hold for pi, 2 = and therefore remain 



valid at least for pi, 2 << 1. This means that, if the anisotropy is weak enough, the leading term in Eq. ^ 
by the dimension A [TV, 0] for any Sn- For all the special cases studied in Sec. this hierarchy persists a 
values of the anisotropy parameters, and the contribution with A [ TV, ] remains the leading one for such TV and pi ,2 



9L]j IS given 
so for finite 



(iii) Of course, it is not impossible that the inequalities (81a), (81b) break down for some values of n, d and pi,2, 
and the leading contribution to Eq. ( pO| ) is determined by a dimension with N' ^ N and/or p > 0. 

Furthermore, it is not impossible that the matrix ( |79| ) for some pi.2 had a pair of complex conjugate eigenvalues, 
A and A*. Then the small mr behavior of the scaling function ^(mr) entering into Eq. (|90| ) would involve oscillating 
terms of the form 



(mr)^'' |Ci cos [im A (mr)] + C'2 sin [im A (mr)] | , 



with some constants Ci. 

Another exotic situation emerges if the matrix (79) cannot be diagonalized and is only reduced to the Jordan form. 
In this case, the corresponding contribution to the scaling function would involve a logarithmic correction to the 
powerlike behavior, (mr)^ [Ci In(TOr) + C2] , where A is the eigenvalue related to the Jordan cell. However, these 
interesting hypothetical possiblities are not actually realized for the special cases studied above in Sec. 0. 

(iv) The inclusion of the mixed correlator (v/) oc nS{t — t')C'{r/€) violates the evenness in n and gives rise to 
nonvanishing odd funct ions 5 2n+i and to the contributions with odd N' to the expansion (po| ) for even functions. If 



the hierarchy relations ( plaj ) , (Bib 
with A[A^, 0]. If the relations Mc 



hold, the leading term for the even functions will still be given by the contribution 
^ hold, the leading term for the odd function iS'2n+i will be given by the dimension 
A[2n, 0] for n < (d + 2)/4 and by A[2n + 1, 1] for n > (d + 2)/4. Note that for the model with an imposed gradient, 
the leading terms for S'2n+i are given by the dimensions A[2n + 1, 1] for all n; see [|4|. This can be related to the 
observation of Ref. |^ that the odd structure functions of the velocity field appear more sensitive to the anisotropy 
than the even functions. 
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Representations similar to Eqs. (|49|), ( |9C| ) can easily be written down for any equal-time pair correlator, provided 
its canonical and critical dimensions are known. In particular, for the operators F[N,p\ in the IR region (Ar — > oo, 
mr fixed) one obtains 

(i^[7Vi,Pi]F[iV2,P2]) = ^ (Ar)-^i"-i'^[«'..'iC^,p;jv',p'(mr), (91) 

N,p N',p' 

where the summation indices N, N' satisfy the inequalities N < Ni, N' < N2, and the indices p, p' take on all 
possible values allowed for given N, N' . The small mr behavior of the scaling functions S,n ,p;N ' .p' {fn^) has the form 

(,N,p;N',p'{mr)= J2 CN"^p"{mr)^^^"'P"\ (92) 

N",p" 

with the restriction N" < N + N' and corresponding values of p"; Cn" ,p" are some numerical coefficients. 

So far, we have discussed the special case of the velocity correlator given by Eqs. (||), (p^). Let us conclude this 



Section with a brief discussion of the general case (|ll| ) . The RG analysis given above in SecsT^I|-[VI| can be extended 
directly to this case; no serious alterations are required. From the expressions (pi]), ( [zi] ) it immediately follows that 
only even polynomials in the expansion ( [l^ ) can give contributions to the renormalization constants, and consequently, 
to the coordinates of the fixed point and the anomalous dimensions. For this reason, the odd polynomials were omitted 
in Eq. (^2|) from the very beginning. Moreover, it is clear from Eq. ( ^ ) that only the coefficients a; with Z = 0, 1 
and bi with Z = 0, 1,2 contribute to the constants Z12 in Eq. ( p8[ ) and therefore to the basic RG functions ( p3[ ) and 
to the coordinates of the fixed point in Eq. (^2[). Therefore, the fixed point in the general model ( [Tl|) is parametrized 
completely by these five coefficients; the higher coefficients enter only via the positivity conditions (|13|). 

Furthermore, it is clear from Eqs. (|7^) that for a = 0, only coefficients a/ with I < 2 and 6; with I < 3 can contribute 
to the integrals i7„ and, consequently, to the one- loop critical dimensions (^). Therefore, the calculation of the latter 
essentially simplifies for the special case oo = 1, ai = and hi — for Z < 2 in Eq. ([ll|). Then the coordinates of 
the fixed point (^2|) are the same as in the isotropic model, in particular, a* = 0, and the anomalous exponents will 
depend on the only two parameters 02 and 63. We have performed a few sample calculations for this situation; the 
results are presented in Figs. 6-9 for A[n,p] with n = 3, 4, 5, 6 vs 02 for 63 = 0, vs 02 — b^, and vs 63 for 02 — 0. In 
all cases studied, the general picture has appeared similar to that outlined above for the case (|l^). In particular, the 
hierarchy of the critical dimensions, expressed by the inequalities (|8^), persists also for this case. We may conclude 
that the special case (p^) case represents nicely all the main features of the general model (11). 



VIII. EXACT SOLUTION FOR THE SECOND-ORDER STRUCTURE FUNCTION AND 

CALCULATION OF THE AMPLITUDES 



In renormalized variables, dimensionality considerations give, see Eq. (^4|) 

S2ir) = iryiy)Rif,r,mr,g,a), (93) 

where R = R2 is some function of dimensionless variables; the dependence on c?, e, pi^2 and the angle between the 
vectors r and n is also implied. The asymptotic behavior of the function ( p3| ) in the IR region (Ar >> 1, mr fixed) 
is found from the solution of the RG equation, see Eqs. (^^ and (|48| ) in Sec. 0: 

S2{r)=D^'r^~'g,a^r), (94) 

where the scaling function (J(u) is related to the function R from Eq. ( p3[ ) as follows: 

^{u) = R{l,u,g^,a^), (95) 

with 5, and a* from Eq. (^). The function R in Eq. ( |9^ ) can be calculated within the renormalized perturbation 
theory as a series in g, 



00 



^(ff,---) = V5"^n(---), (96) 



n=0 



where the dots stand for the arguments of R other than g. Making the substitutions fir ^ 1, g ^ g^t and a a^,, 
expanding i?„ in e, and grouping together contributions of the same order, from Eq. ( |9^ ) we obtain the s expansion 
for the scaling function: 
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au)^y^e^^Uu)- (97) 



n=0 



[It is important here that the calculation of any finite order in e requires only a finite number of terms in Eq. ( pM 
because 5, = 0{e), and the coefficients i?„ for the renormalized quantity do not contain poles in e.] However, the 
expansion ( |97| ) is not suitable for the analysis of the small u behavior of ^ (it) : we shall see below that the coefficients 
contain IR singularities of the form uPItl^u, these "large IR logarithms" compensate for the smallness of e, and 
the actual expansion parameter appears to be e In m rather than e itself. 

The formal statement of the problem is to sum up the expansion ( p^ ) assuming that e is small with the additional 
condition that e\nu — 0{1). The desired solution is given by the OPE; it shows that the small u behavior of £^{u) for 
the second-order structure function has the form, see Eqs. (89) and ( |90| ) in Sec. 
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The first term is the contribution of the operators 1 and F[2, 0] with the dimension A[2, 0] — 0, see the text below 
Eq. (^), and A [2, 2] is given to the order 0{e) in Eq. (|8^); the dots in Eq. ( p8|) stand for the contributions of order 
and higher. The coefficients Ci depend on e, d, pi^2, and the angle 1}. They are sought in the form of series 

in e: 

00 

Cn = J2ci'^e\ (99) 

fe=0 

Equation ( |9^ ) can be expanded in e and compared with (p7|). Such a comparison with known e expansions for the 
coefficients ^„ in ( |9^ ) and exponents in ( p8| ) allows one to verify the representation ( p8| ) directly from the perturbation 
theory, and to determine the coefficients in Eq. (^9[). 

In order to find the coefficients Ci,2, related to the exponents of order 0(e) in Eq. (p8|), we do not need the entire 
functions ^n(w), we only need the terms of the form In'' u in their small u asymptotic expansions. Knowledge of the 
terms 1 and e In u is sufficient for the calculation of the lowest (zero) order coefficients C^^g • Calculation of the terms 

(e In u)P with p > 1 gives additional equations for the same quantities c[^l , which must be satisfied automatically 
and can therefore be used to verify the representation (|9^). Calculation of the terms e^'ln'' u with p > q gives the 
equations for the higher terms with k > 1. Finally, the terms of the form e^tt'^ln'u with fc > are related to 
the contributions of the form in Eq. (p8|). 

In what follows, we confine ourselves with the calculation of the coefficients C^fl . To this end, we have to find the 
fimction ^(u) with the accuracy of 

^{u) = Ao + Aielnu, (100) 

where Ai — Ai{d, pi,2,'&)- ^From the definition of the structure functions (||) and the representations (p^), ( p5[ ) it 
follows that ^(u) = {D{uj, k))RG, where by the brackets we have denoted the following operation 

{Dicu, k))nG f^Jj^ Die., k) \l - exp(ikr)l (101) 

and -D(w, fc) — {90) is the exact pair correlator of the field 0, which we need to know here up to the one-loop 
approximation: 

D{uj,k)= + . --'^"r + + r'""-- . + ... (102) 



In what follows, we set ^ = 00 in Eq. (g); the IR regularization is provided by the "mass" m in the correlator 
(^. In the momentum-frequency representation, the correlator (H) takes on the form C(fc) — (27r)''i5(k), and the 
bare propagator in the renormalized perturbation theory is obtained from {69) q in Eq. (Hq) with the substitution 
i^ofc^ e(fc) = vk'^ + ;/Q;(kn)^. 

Therefore, the contribution of the first (loopless) diagram in Eq. (|102|) to the function (|101|) has the form 



1 

2d 



a*(l + 2cos2 1)) 
{d + 2) 



0(P?,2) 



(103) 
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For the sake of simplicity, here and below we confine ourselves with the first order in the anisotropy parameters pi.2 
(we recall that a* = 0{pi^2))', only the final results will be given for general p. Note also that all the integrals entering 
into the left hand side of Eq. (103) are well-defined for £ — oo, for example: 



fc2 



1 — cxp(ikr) 



(5(k) kikj 



1 — exp(ikr) 



2d{d + 2) 



(104) 



It is also not too difficult to take into account the contributions from the second and third (asymmetric) diag ram s 
in Eq. (102). It follows from Eq. ( ^4[ ) and (|35| ) that one simply has to replace e(fc) in the left hand side of Eq. (103) 
with the expression vok'^ + aoi^o(nk)2 — ^o'e and change to the renormalized variables using the relations (29) and 
(^). The pole part of the integral J{m) in Eq. (|3^) is then subtracted; the e expansion of the difference gives rise to 
the logarithm of m. The substitution g ^ g^,, a ^ a* leads to drastic simplification (this is due to the fact that the 
expressions (^ ) for th e fixed point originated from the sa me E q. (|35|)), and the total contribution of the first three 
diagrams in Eq. (102) reduces to the loopless cont ribution ( |103| ) multiplied with the factor (1 — elnu). 

Now let us turn to the last diagram in Eq. (102). The quantity to be "averaged" in Eq. ( |lOl[ ) has the form 



duj' f dq gvp!^Tij{c() 



C{q)qtqj 



where q' = k — q. We perform the integrations over uj' and lo (the latter is implicit in Eq. (101)) and replace q' with 
k in the integral over q (we recall that C(fc) cx ^(k)); this gives 



4e(fc)(fc2 + m2)'i/2+'^/2 



(105) 



Using Eq. (O) and the relation 



dqS{q) 



i/d 



vd{d + 2) 



2ninj 



we obtain for the quantity ( |105| 



Advk"^ (fc2 + m2)'i/2+'^/2 



with the coefficients 

CO = 1 + P2 - a(<i+ l)/(fi + 2), ci 



Co + Ci cos ip 



{d - l)pi ~ p2~ a{d^ + d- 4)/(d + 2) 



(106) 



(we also recall that ip is the angle between k and n). Substituting Eq. (105) into (101) gives 

|^[co/o(m) +ci/i(u)] 

with the dimensionless convergent integrals 



In(u) = -J 



dk 



r2 J (277)^^ fc2(fc2 +m2)'i/2 



1 — exp(ikr) 



(107) 



(108) 



[within our accuracy, we have set e = in all exponents in Eqs. (107), (108)]. The leading terms of their small u 
behavior are found as follows: one differentiates /„(m) with respect to u, performs the change of variable k uk, 
expands in u the quantity in the square brackets, and selects the leading terms of order 1/u. This gives 



Io{u) - Inw, hiu) 



Cd In?! 



2d(d + 2) 



1 + 2 cos^ d 



(109) 



with Cd from Eq. 



Collecting all the contributions of the form (100) in Eq. (101) within our accuracy gives 

a,(l + 2cos2i9)1 {d+\)[{d'^+d-A)pi~dp2] 



^0 



1 

2d 



d+2 



(d- l)2(d + 2)3 



P2(C0S?9), 



(110) 
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where P2{z) = — 1 / d is the second Gegenbauer polynomial. 

It is sufHcient to take A = A[2, 2] in the isotropic approximation, A — ed{d + \)/{d— l){d + 2), see Eq. (|8^) [it is 
important here that Ax — O(/0i,2)L 

Comparing Eq. (98) with Eq. (IOC) gives C2 = Aie/A and Ci = — C2. Substituting these expressions into Eq. 
(pS) and then into Eq. (B4) and using Eq. we obtain the final expression for 82- 



S2{r) 



{d-l)Cd 



1 - (p, + p2)/d + 2%::# P2(cos ^) - {mrr + ^ ' /Vi -/P^] ^^^^^^ 



(d + 2) 



(d-l)(d + 2)2 



(111) 



Similar (but much more cumbersome) calculations can also be performed if the anisotropy parameters pi^2 are not 
supposed to be small. The result has the form: 



^2(r) = 



deD;;^r^-'' 



-(mr) 



{d - l){d + pi + p2)Cd 
^ (d + a^) 



1 



.g*(d + a*)(rf+ 1)(P2 - Pi) 



F2(cosi9) 



{d-l)d{d + 2) 



P2(cos?9) 



a^SdF* - {d + 2)F* 



2{d- l)(d + 2)A 

g,{d+l){p2- pi) 



2A 



?>{d + a,) 
~ {d + 2) 

3{d + a*) 
(d + 2) 



F* 



(112) 



where the quantities 5,, a, and A = A[2,2] are given in Eqs. ( ^2[ ) and (|85|), ^2(2) = — (1 + a,)/((i + a,) and 
F* = F(l, 1/2 + n; d/2 + n; -a,) (as in Eq. (H)). 

Now let us turn to the exact equation satisfied by the second-order structure function S'2(r). /^From the definition 
of the latter and Eq. (E3) we obtain 



2iyoAS2(r) + DoS^,{r)d,djS2{r) = 2C{r/e). 
In the inertial range we neglect the viscosity and set C{r/£) ~ C(0) — 1; this gives: 

S,,{r)d,djS2{r)^2D^\ 



(113) 



(114) 



We should also set mr << 1. The integral (|2^) has a finite limit for m = 0; the explicit calculation with Tij from 
Eq. (|l|) gives: 



2r^ 

in - pi) 



{d + e — l)5ij — eriTj/r"^ + pi(l + ez^) + p2 {d + e — 2)ninj — ez{nirj + njri)/r 



6ij + 2ninj + e{5ijz'^ + Tirj/r^ + 2z{jiirj + njri)/r) + e{e — 2)z'^rirj /r^ 



id + e + 2) 
where z = (nr) /r and 

C{e) = -2^+^ (4^)'^/2 r(rf/2 + e/2 + l)/r(-e/2) = 2ed/Cd + 0{e^). 



(115) 



(116) 



For what follows, we need to know the action of the differential operator in the left hand side of Eq. (114) onto a 
function of the form r'^~'^~^'^ip{z), where 7 is an arbitrary exponent (below we use 7 to denote different exponents; the 
precise meaning in each case will be clear from the context): 



^Piz) 



C{e) 



/:(7)^(z) 



(117) 



Here £(7) is a second-order differential operator with respect to z, whose explicit form is 

£(7) = (2 7 - e){(rf - l)(d + l){d + 2 + e) + {d- l){d + 2)(pi + P2) + 2e{dp2 - Pi) + z^ed{d + l)(pi - P2) + 
+ ^[{d+ l)pi + (e + l)p2 + z^{pi{ed -e-2))+ p2{d^ - 2 - e)] } + |(1 - d){d - 1 + e){d + 2 + e)- 

- {d^ + l)pi + {d^-d^ \)p2 - e{d - l)pi - eid^ - 2)p2 + 4eVi - (2d + l)e^P2 + z^ [2pi - [d^ - 2)p2 + 
+ e{d + 1) ((1 - d)pi + {2d - 3)p2) + {-{d + 3)pi + {2d + l)p2)] - 7(1 - ^') [4pi - 2{d^ - 2)p2 + 
+ e(4pi - 2dp2)] I 2(9^ + |(d + e - l)[d + 2 + e + pi + (d+ 1 + e)p2] + [(2 + e(d + 1 + £))pi - 



{{d + . 



2-e 



)P2]}{1-Z^)d' 



(118) 
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We seek the solution to Eq. (114) in the form 

which gives the following differential equation for ')jj{z): 

Cip{z) = l, £ = £(7 = 0). 

Here the operator C has the form 

£ = e2{z)d^,+ii{z)d, + eQ(z), 



(119) 
(120) 
(121) 



with some functions £i{z) whose explicit form is readily obtained from Eq. ( 118 ) and will not be given here for the 
sake of brevity. The operator C is self-adjoint on the interval — 1 < 2: < 1 with respect to the scalar product 



with the weight function 



The integration gives 



if, 9) = I ^dzp{z)f{z)g{z) 



p{z) = £^^{z) exp j dz'h[z')ll2{z'). 



n(.\ - n - .2Ud-3)/2 , [2 + £(d + 1 + £)]pi - {{d + g)^ - 2 - e\p2 , 

^^ ^"^ ^ \ (d + e-l)[2 + d + e + pi + (d+l + e)p2] 



where 



(2 + - £)[-(l + £)pi + {d{d + £)/2) - l)p2] 

[2 + e(d + 1 + e)\pi + [2 + £ - (rf + eY\p2 



(122) 



(123) 



(124a) 



(124b) 



(the lower limit in Eq. (123) is chosen such that p{z = 0) = 1). 

The self-adjoint operator £ possesses a complete set of eigcnfunctions, but it is not possible to find them explicitly 
in the general case. However, Eq. ( |119| ) can be treated within regular perturbation expansions in e, p\^2 and Xjd. In 
what follows, we shall discuss them separately. 

Perturbation theory in e. We write the operator £ — £(7 — 0) in Eq. (IIS) and the desired solution '4>(z) as series 
in e: 



£ = £0 + £ £e 



^{z) = V^°^(^)+£V'^'^^) + 



The leading operator £0 has the form (121) with the functions 

l,{z) = I,{z) \{d - \){d + 2) + pi(d + 1) + P2] , 

where 

7o = 2(d + a»), =-z[(i-l-a,(l-z2)], Z2 = (1 - z^)[l + a,(l - z^)] . 

It is easy to see that the constant is one of the eigcnfunctions of the operator £0 (a consequence of the fact that Iq 
is independent of z). The other eigenfunction (with the eigenvalue of zero) is P2{z) from (112): 

£oP2(^) = 0. 

Therefore, a solution to the equation Caip^^'^iz) = 1 can be written in the form 

^'^''\z)=eo^ + CP2{z). (125) 
In order to find the constant C, consider the equation for the correction -0'^^^ (z): 
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(126) 



Using the scalar product (122) with the weight (|l2J) at e = 0, 

p(0)(2) = (l-z2)(<i-3)/2 



(127) 



(the operator Co is self-adjoint with respect to it) we obtain {P2, Cq^''-^^) — (£o^^j V^'-"'^-') = 0- Then from (126) it 
foUows that (P2,'Ce'i/'*"^) = 0. From the latter relation and Eq. ( p5| ) we find 



c = -- 



which gives 



*0 



(128) 



Below we shall show that the sol ution (128) coincides with the RG result (112) for m = 0. To this end, we shall 
show first th at th e last term in Eq. (112), proportional to (mr)^, is an eigenfunction of the operator in the left hand 
side of Eq. (113) with the eigenvalue of zero. We seek the z ero mode of that operator in the form 7-2-e+7 $(2) with 
some exponent 7 and function $(z). According to Eq. (|ll7|), we arrive at the equation £(7)$(z) = with £(7) from 



( |118| ). For e = 0, the solution of the latter equation is given by 7 = and <i>(z) = P2iz) from Eq. (|112| ). The 0(e) 
correction to the solution 7 = is sought from the equation 

{£0 + eC, + jC^)^z) ^ 0, 



where jC^ is the part of operator (118) linear in 7. We require that the correction operator eC^ + ^C-y do not shift 
the original (zero) eigenvalue. This leads to the equation {P2, {eC^ + ^Cy)P2) = 0, from which it follows 

,{P2,C,P2) 



7 = -e 



iP2,C^P2) 



(129) 



Using Eqs. (118) and (127) it is not difficult to check that the expr ession (129) coincides with the dimension (|85|), 
e.g., 7 = A[2, 2]. Therefore, we have shown that the last term in (112) is indeed an eigenfunction of the operator from 
(114), and the exponent and function ^{z) are found in the first and zeroth order in e, respectively. 
Using Eq. (|129|), expression (|128D can be rewritten in the form 



e{P2,C,-l)^ 
A[2,2] {P2,C-fP2) 



P2{Z) 



(130) 



Calc ulation of the scalar products using Eqs. (118) and (|127| ) shows that Eq. (|130| ) indeed coincides with the solution 
( pl^ ) for m = 0. 

Perturbation theory in the anisotropy parameters pi,2- In this case it follows from Eq. (118) with pi^2 = 0, that the 
leading operator £0(7) has the form (121) with 

4 = (2 + 7-e)(d- l)(d + 7)(d + 2 + e), 

11 = -{d - l){d - 1 + e){d + 2 + e)z, 

12 = ((i-l+e)(d + 2 + £)(l-z2). 



(131) 



It can be represented in the form 



£0(7) = -(d - 1 + £)((i + 2 + £)£ + fo- 



(132) 



The eigenfunctions of the operator C 



1)9^ + {d— l)zdz are nothing other than the well known Gegenbauer 



polynomials P„(z), the corresponding eigenvalues being n(n + d — 2); see pT} - They form an orthogonal sys tem on 
the interval [-1, 1] with the weight (1 - z2)('^-3)/2 jg^^ (p4|) with pi,2 = 0]. Since the term £0 in Eq. (|l3|) is 
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independent of z, the eigenfunctions of the operator £0(7) are the same polynomials P„(z), and the eigenvalues have 
the form: 



{d-l + e){d + 2 + s)n{n + d - 2) + ig, 



(133) 



with io from Eq. ( |131| ). 

The knowledge of the eigenfunctions and eigenvalues of the operator Co = £0(7 = 0) in Eq. ( |120D allows one to 
easily construct the iterative solution of that equation. In particular, in the leading order one has 

^P^°\z) = 1/io = l/(2-e)d(d- l)(d + e). 



For TO 7^ 0, the corrections to the solution of Eq. (114) are determined by the zero modes of the operator £(7). In 
the leading (zeroth) order in pi,2 they are simply found by equating the eigenvalues ( |133| ) to zero. The "admissible" 
(in the sense of l^^) solution is: 



27,1"^ = {e-d-2) + y/{e + d- 2)2 + 4n{n + d - 2)(d + e - l)/{d - 1). 



(134) 



Therefore, in such approximation the zero modes of the operator from Eq. (114) have the form r^^^+'^'n 'p„(z), in 
agreement with Ref. (owing to the z ^ —z symmetry of our model, only even n can contribute). From the RG 
viewpoints, the exponents in Eq. ( |134| ) are related to the 71-th rank composite operators built of two fields 9 and n 
derivatives; in particular, 72 (which gives the leading correction for mr << 1) coincides with the dimension A[2, 2]. 

In order to find the 0(pi,2) correction to 72"'', which we denote 72^"*, we expand the operator £(7) from Eq. (11?) 

around the point pi 2 = and 7 — j^"^ up to the terms linear in pi 2 and 72^'': 



£(7)^£o + 5]p.r + 7^"A, 



which defines the operators £* and £1. Proceeding as for the previous case, we obtain 



(-P2,£lP2)' 



where the scalar product is defined with the weight (1 — z^)*^"^ 3)/'^^ The explicit calculation gives the desired result: 



72 



(1) 



4{d-2){d+l)e [2pi + (rf + £)p2] 



(d - l)(d + 4)(d + 2 + e) [4d(d2 - i)e +{d- l f{d + 2 - e)2]i/2 ' 
Perturbation theory in 1/d. For d — > 00, the weight ( |l24[) takes on the form 

d/2 



p(z) 



1 



1-P2ZV(1+P2) 



(135) 



(136) 



^From Eq. (13(:) it follows that, for d » 1, the function p(z) differs remarkably from zero only in a small vicinity 
of the point z ~ 0, more precisely, z^ < i/^f Th is becomes especially clear if one changes to the new variable 



(137) 



X = dz /2(1 + P2)- In terms of the fmiction ( |136| ) has a finite limit for d 



The integration domain in the s cala r product (122) becomes — -v/d < x < \/d, which for d ~> 00 gives —00 < x < 00. 

According to ( |118D , equation (12C) in terms of the variable x and in the leading approximation in 1/d can be written 
in the form 



where 



Coi^ix) = 2/d\ 
Co = £0(7 = 0), £0(7) ^dl- 2xd, + 2(2 + 7 - e). 



(138) 
(139) 
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The eigenfunctions of the operator — 2xdx are the well-known Hermit polynomials Hn{x) (see, e.g. 

{dl~2xd,)H„{x) ^2nHn{x), 71 = 0, 1, 2, • • • , 



which form an orthogonal system on the axis —00 < a; < 00 with the weight (137). It is clear from Eq. (139) 
that the eigenfunctions of the operator £0(7) are also the polynomials iJ„(a;), the corresponding eigenvalues being 
A„ = 2(n + e-2-7). 

The zero modes (A„ ~ 0) correspond to the values 

+ (140) 

where n — 2, A, 6 and so on, owing to the symmetry z — > ~z of our model; the minimal e igen value is 7 = e with 
n = 2. Therefore, the operator Cq = >Co(7 = 0) has no zero modes, and the solution to Eq. (13S) is simply given by 
V'(x) = l/d3(2-e). 

The corrections to the solution 7 = £ in 1/d are found as above for the perturbation theories in e and pi^2- This 
gives: 

7 = £ + 2£(1 + 2p2 + PlP2)/(f + 0{l/d^) 

[the 0{\/d) term is absent], and the corresponding zero mode is 

(1 - x") - p2{2 + PI+ p2)/d{l + P2) + 0{lld^). 



IX. CONCLUSION 



Let us recall briefly the main points of the paper. 

We have studied the anomalous scaling behavior of a passive scalar advected by the time-decorrelated strongly 
anisotropic Gaussian velocity field. The statistics of the latter is given by Eqs. (||) and (ITTl). The general case ( |ll| ) 
involves infinitely many parameters; most practical calculations have been performed for the truncated two-parametric 
model (14), which seems to represent nicely the main features of the general case (|lj). 

The original stochastic problem (0)-(|^) can be cast as a renormalizable field theoretic model [Eqs. (p^)-(p^], which 
allows one to apply the RG and OPE techniques to it. The corresponding RG equations have an IR stable fixed point, 
Eq. dH), which leads to the asymptotic expressions of the type (^, ( ^ for various correlation functions in the 
region Ar >> 1. 

Those expressions involve certain scaling functions of the variable mr, whose behavior at mr << 1 is determined 
by the OPE. It establishes the existence of the inertial-range anomalous scaling behavior. The structure functions are 
given by superpositions of power laws with nonuniversal (dependent on the anisotropy parameters /Ci,2) exponents. 

The exponents are determined by the critical dimensions of composite operators (pq) built of the scalar gradients. 
In contrast with the isotropic velocity field, these operators in our model mix in renormalization such that the matrices 
of their critical dimensions are neither diagonal nor triangular. These matrices are calculated explicitly to the order 
0(e) [Eqs. ( [t^, ( [76|), ([78|), ([79|)], but their eigenvalues (anomalous exponents) can be found explicitly only as series 
in pi^2 [Eqs. (^6[), (|87[)] or numerically [Figs. 1-9]. 

In the limit of vanishing anisotropy, the exponents can be associated with definite tensor composite operators built 
of the scalar gradients, and exhibit a kind of hierarchy related to the degree of anisotropy: the less is the rank, the 
less is the dimension and, consequently, the more important is the contribution to the inertial-range behavior [see 
Eqs. (|l|)].^ 

The leading terms of the even (odd) structure functions are given by the scalar (vector) operators. For the finite 
anisotropy, the exponents cannot be associated with individual operators (which are essentially "mixed" in renormal- 
ization), but, surprising enough, the aforementioned hierarchy survives for all the cases studied, as is shown in Figs. 
2-9. 

The second-order structure function S'2(r) is studied in more detail using the RG and zero- mode techniques; like 
in the isotropic case its leading term has the form ^2 cx: r^~^, but the amplitude now depends on pi^2 and the 
angle between the vectors r and n from Eq. (p^. The first anisotropic correction has t he fo rm (mr)^!^'^^ with 
the exponent A[2, 2] = 0{e) from Eq. (|85|). The function ^2 satisfies the exact equation ( |113| ); this allows for an 
alternative derivation of the perturbation theory in £ , and also gives a basis for perturbation theories in 1/d or the 



anisotropy parameters pi^2] see the discussion in Sec. VIII 



It is well known that, for the isotropic velocity field, the anisotropy introduced at large scales by the external forcing 
or imposed mean gradient, persists in the inertial range and reveals itself in odd correlation functions: the skewness 
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factor Ss/Sl^^ decreases for mr but slowly (see Refs. PP,p|,p^), while the higher-order ratios 5'2„+i/5'2^^^^ 
increase (see, e.g., p^ , [l6| , p^ ) . 

In the case at hand, the inertial-range behavior of the skewness is given by Ss/S^^^ cx (mr)^[3'il. For pi,2 0, 
the exponent A [3, 1] is given by Eq. (jsy) with n — 3 and p = 1; it is positive and coincides with the result of Ref. 
1^. The levels of the dimension A[3, 1] on the pi-p2 plane are shown in Fig. 10. One can see that, if the anisotropy 
becomes strong enough, A[3, 1] becomes negative and the skewness factor increases going down towards to the depth 
of the inertial range; the higher-order odd ratios increase already when the anisotropy is weak. 
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TABLE I. Canonical dimensions of the fields and parameters in the model (|26[ 
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FIG. 1. Levels of the dimension A[2, 2] for d = 3 on the plane pi-p2- Value changes from 1.15 (left-bottom) to 1.4 (right-top) 
with step 0.05. 




FIG. 2. Behavior of the critical dimension A[3,p] for d = 3 with p = 1,3 (from below to above) vs pi for p2 = — left, vs 
p = pi = p2 — center,vs p2 for pi = — right. 




FIG. 3. Behavior of the critical dimension A[4,p] for d = 3 with p = 0, 2,4 (from below to above) vs pi for p2 = — left, vs 
p = pi = p2 — center, vs p2 for pi = — right. 




FIG. 4. Behavior of the critical dimension A[5,p] for d = 3 with p = 1, 3, 5 (from below to above) vs pi for p2 = — left, vs 
p = pi = p2 — center,vs p2 for pi = — right. 
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FIG. 5. Behavior of the critical dimension A[6,p] for d = 3 with p = 0, 2, 4, 6 (from below to above) vs pi for p2 = — left, 
vs p = pi = p2 — cx.nter,vs p2 for pi = — right. 




FIG. 6. Behavior of the critical dimension A[3,p] for d = 3 with p = 1, 3 (from below to above) vs a2 for 63 = — left, vs 
02 = 63 — center,vs &3 for a2 = — right. 




FIG. 7. Behavior of the critical dimension A[4,p] for d = 3 with p = 0, 2,4 (from below to above) vs a2 for 63 = — left, vs 
0.2 = bs — center,vs 63 for a2 = — right. 




FIG. 8. Behavior of the critical dimension A[5,p] for d = 3 with p = 1, 3, 5 (from below to above) vs 02 for 63 = — left, vs 
(12 = bs — center,vs 63 for 02 = — right. 




FIG. 9. Behavior of the critical dimension A[6,p] for d = 3 with p = 0, 2, 4 (from below to above) vs a2 for 63 = — left, vs 
0.2 = bs — center,vs bs for 02 = — right. 
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